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ABSTRACT 

HD 181068 is the brighter of the two known triply eclipsing hierarchical triple stars 
in the Kepler field. It has been continuously observed for more than 2 years with the 
Kepler space telescope. Of the nine quarters of the data, three have been obtained in 
short-cadence mode, that is one point per 58.9 s. Here we analyse this unique dataset 
to determine absolute physical parameters (most importantly the masses and radii) 
and full orbital configuration using a sophisticated novel approach. We measure eclipse 
timing variations (ETVs) , which are then combined with the single-lined radial velocity 
measurements to yield masses in a manner equivalent to double-lined spectroscopic 
binaries. We have also developed a new light curve synthesis code that is used to 
model the triple, mutual eclipses and the effects of the changing tidal field on the stellar 
surface and the relativistic Doppler-beaming. By combining the stellar masses from the 
ETV study with the simultaneous light curve analysis we determine the absolute radii 
of the three stars. Our results indicate that the close and the wide subsystems revolve 
in almost exactly coplanar and prograde orbits. The newly determined parameters 
draw a consistent picture of the system with such details that have been beyond reach 
before. 
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1 INTRODUCTION 



The 



Kepler space telescope, in addition to its primary 
science aims, has led to a new era in the investiga- 
tion of multiple star systems. Among the highlights we 
find t he discoveries of th e first triply eclipsing triple sys- 
tems jCarter et al.l I2OIII : [Dcrckas ct al. 201 1j) and some 
inter e sting studies of multiple star syst e ms Ipteffen et al] 
2OIII: iFeiden. Chabover. fc Dotted I2OIII : iGies et al.l |2012| : 
Lehmann et al.ll2012l '). 

Binary and multiple systems have an important role 
in astrophysics. The most accurate way to measure stellar 
parameters is through eclipsing binaries, and their distance 
determination is also very accurate. Their light curves pro- 
vide essential information on the internal structure of the 
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components, their atmospheres and their magnetic activity. 
In the case of noncircular orbits and multiple systems, the 
orbital elements can change significantly, allowing detailed 
insight into the time variation of these parameters. The spe- 
cial geometry of the very rare and new category of eclipsing 
systems, namely the triply (or mutually) eclipsing triple sys- 
tems, enables us fast and easy determination of further char- 
acteristics that otherwise could only be studied with great 
effort on a long time-scale. 

As an example, we refer to the spatial configura- 
tion of such hierarchical triple systems, which is a key- 
para meter in understanding the ir origin and evolution (see 
e. g. iFabrvckv Sz Tremaind 20071 . and references therein). In 
the absence of mutual eclipses, the two ways to determine 
the mutual (or relative) inclination in a hierarchical system 
are (a) astrometric (or, more rarely, polarimetric) measure- 
ments of the spatial orientations of the two orbits individu- 
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ally, or (6) indirect dynamical calculation from the measured 
mutual gravitational perturbations of the bodies. The first 
method requires long baseline optical (or very-long baseline 
radio) interferometric measurements for the most interest- 
ing close binaries, which typically have milli-arcsecond an- 
gular separations. It is there fore not suprising tha t, starting 
with the pioneering work bv lLestrade et ahl (|l993l ) on Algol, 



naries (see also iBaron et al.ll2012l: 


Sanborn & Zavalall2012l: 


IPeterson et al.ll2011; O'Brien et all 


2011. for more recent re- 



suits). The applicability of polarimetric measurements (al- 
though does not require high-catego ry instrumen ts) in this 
field is even more restricted (see e. g. |Piirolall2010l ). The sec- 
ond method, the detection of gravitational perturbations, re- 
quires accurate, frequent and continuous photometric eclipse 
time determination. This method will be described in detail 
in the next section. 

The situation is much easier in the case of mutual 
eclipses, where the shape of the light curve (especially 
around the ingress and egress phases) contains direct and 
unique information about the system geom e try. T his is dis- 
cussed in detail by iRagozzine fc HolmanI (|2010l ) and IPall 
l|2012l ). The former authors list several other values of 
multi-transiting systems, mainly in the context of multi- 
ple planetary systems. Their model has been succesfuUy 
applied to analysing complex light curves and determin- 
ing the corresponding geometrical and physical parameters 
(both for the orbits and the i ndividual bodies) for different 

Kepler- 



mu ltiple-transiting p lanetary (iLissauer et al. 



2011 



ll. lDovle et al.ll201ll. Keple r-16, IWelsh et al.ll2012l . Kepler- 
34b-35b. ICarter et al.|[20l3 . Kepler-36) and stellar systems 
l|Carter et al.ll201ll . KOI-126). 

KOI-126 and HD 181068 are the first representatives of 
the new category of the triply eclipsing triple systems. Both 
are also members of a very small group of compact hier- 
archical triple stellar systems. They contain a close binary, 
with orbital periods °^ = 1.77 and = 0.91 days, and 
a more distant component forming a wider binary with the 
centre of mass of the close pair with periods PI'^' = 33.92, 
and P™ — 45.47 days, respectively. The main speciality of 
the two systems is their triply eclipsing nature, which means 
that both the inner and the outer binaries show eclipses. 
They have other, very peculiar characteristics. Both belong 
to the most compact triple stellar systems, and there is only 
one known hierarchical triple system with a shorter outer 
period, namely A Tan, with P2 = 33.03 days. Furthermore, 
these two systems are unusual even amongst the very few 
similarly compact triples, in having reversed outer mass- 
ratio. In other words, in these two objects the wide, single 
component is the more massive star, and also the largest and 
brightest. Before Kepler, the highest known outer mass-ratio 
did not reach 1.5, and for 97% of known hierarchical triplets 
it remained under 1, i. e. almost in all the catalogized sys- 
tems, the total mass of the clo se binary exceed ed the mass 
of the tertiary component (see lTokovinin|[2008l ) . (The ques- 
tion of whether this comes from observational bias is not 
discussed here.) In contrast, the outer mass ratios of these 
two new systems aregAe' ~ 3-0, and g™ ~ 1-9, respectively. 

Despite the similarities of KOI-126 and HD 181068 to 
each other, there are remarkable differences between the two 
systems. On one hand, KOI-126 consists of three nearly 
spherical main sequence stars, where the members of the 



close binary have such a low surface brightnesses that their 
light curve modelling is largely equivalent to those of the 
multiple planetary systems. This is not true for HD 181068, 
where all the three stars are tidally distorted, have almost 
equal surface brightnesses and show evidence of intrinsic 
light variations, all of which make light curve modelling of 
HD 181068 more difficult than for KOI-126. On the other 
hand, dynamical analysis of HD 181068 is much less complex 
than for KOI-126, because of the much simpler and appar- 
ently constant orbital configurations. As a consequence, our 
method of light-curve analysis is much closer to the tradi- 
tion al eclipsing binary star light curve modeling methods 
fsee iKallrath fc Milone|[2009l . for a review) than the proce- 
dures applied for systems like KOI-126. 

In this paper, we analyse more than 2 years of Kepler 
observations of HD 181068. We mainly concentrate on de- 
termining the fundamental astrophysical parameters of the 
three stars and orbital elements of the close and wide orbits. 
These quantities by themselves carry very important infor- 
mation already about the system and their members' origin 
and evolution and, furthermore, give the necessary input 
parameters for other forthcoming studies, for example for a 
comprehensive study of pulsations of the red giant compo- 
nent. Nevertheless, due to the uniqueness of the studied sys- 
tem, our aim is not simply to give a case study. The specifics 
of HD 181068 allow us to present methods never used be- 
fore. For example, in our period study (Section 3) , which de- 
pends on the analysis of the eclipse timing variations (ETV) 
for both the close and the wide systems, we determine the 
(inclination-dependent) masses of the wide binary members 
in a new manner. While the radial velocity curve of the most 
massive A component is known, the missing second radial 
velocity curve of the spectroscopically unseen B component 
(i. e. the close binary itself) is replaced by the light-time or- 
bit of the B component deduced from the ETV analysis of 
the shallow eclipses. This method is fundamen tally different 
from the one followed bv ISteffen eFal] (I2OIII ) for KOI-928, 
for example, because it does not use the dynamical part of 
the ETV, only the simple geometrical light-time contribu- 
tion. More details are given in Section 3. In Section |4j the 
light curve analysis procedure is described in detail, while 
Section 5 contains the discussion of the results. Finally, the 
details of our light curve synthesis and analysis code, and 
some additional examples of calculations of certain quanti- 
ties purely in a photometrical and geometrical way from the 
mutual eclipses, are given in the appendices. 

It is i mportant to establis h a clear notation for this 
system. In lDerekas et al.l (|201ll ) the three components were 
labelled A, B, and C (in order of decreasing masses and lu- 
minosities). Here, we use the more clarified and expressive 
denotations, A, Ba, Bb. As before, A denotes the most mas- 
sive and luminous component (the main component of the 
wider A — B binary) , while Ba and Bb refer to the members 
of the close binary formed by the two red dwarfs, formerly 
denoted by B and C. When referring to any physical quanti- 
ties of the individual stars, we use subscripts. For example, 
rriA and mea denote the masses of the A and Ba compo- 
nents, respectively, but me refers to the total mass of the 
close binary, i.e. (mBa + ?TiBb), and rriAB stands for the total 
mass of the hierarchical triple. With this notation we can 
avoid the confusion with the indices of the orbital param- 
eters of different orbits used for the period study. Namely, 
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following the common usage, the elements of relative orbit of 
the Bb component around its companion, Ba is subscripted 
with 1, whereas the relative orbit of the ternary compo- 
nent A, around the center of mass of the Ba — Bb subsys- 
tem (symbolically represented with B) is associated with the 
subscript 2. However, in terms of light-time and the radial 
velocity, the absolute orbit (i.e. the orbit of some star around 
the center of mass) is to be considered, rather than the rel- 
ative orbits. In these cases, those absolute orbital elements, 
which numerically differ from the corresponding relative or- 
bital element, were naturally denoted by the alphabetic sign 
of the given star, or subsystem. 



2 OBSERVATION AND DATA REDUCTION 

The analysis is based on photometry from the Kepler 
space telescope ( Borucki et al.l 20101: jGiUiland etHI l2010l : 
iKoch et al.ll2010l : I Jenkins et aLboiOal lbh. The dataset is 775 
days long, observed in 6 quarters (Q1-Q6) at long cadence 
(time resolution of 29.4 min) and 3 quarters (Q7-Q9) at 
short cadence (time resolution of 58.9 sec). Since HD 181068 
is a ~ 7 magnitude star, it is heavily saturated, resulting in 
charge bleeding. Therefore, the short-cadence observations 
were obtained using a Custom Made Aperture Mask. This 
was uploaded directly to the spacecraft lookup table and 
shaped precisely to match the shape of the target on the 
detector including the bleeding area. 



2.1 Measuring the times of minima 

The 2.1 year-long observations cover ~ 885 orbital cycles of 
the close pair and 17 revolutions of the wide system. Approx- 
imately 10% of the eclipses of the close binary (hereafter we 
refer to them as shallow minima) occur during the eclipse 
events of the wide system (hereafter deep minima) , and can- 
not be observed. Additionally, a few hundred events escaped 
observation due to data gaps. In all, 1177 of the 1770 shal- 
low minima were analysed. The an alysis of these minima was 
quite a complex task. As shown bv lDerekas et al.) l|201ll ). the 
red giant component shows oscillations on a time scale simi- 
lar to the half of the orbital period of the short period binary. 
In addition, there are long term variations, discussed in Sect. 
|4l which slightly distort the shape of the shallow minima, as 
shown in Fig. [T] This distortion has a significant effect on 
the measurement of the exact times of minima. 

To correct for these distortions, we applied the follow- 
ing method in determining the times of minima. We took 
the ± 0.225 days interval around each minimum and fitted 
low-order (4-6) polynomials outside the eclipses. Then we 
corrected each subset, which resulted in a detrended light 
curve. Finally, to determine the times of minima, we fitted 
low order (5-6) polynomials to the lowest parts of the min- 
ima. 

We also analysed the available deep minima. Out of the 
34 events we were able to determine times of minima in 28 
cases. (One of these events was omitted from the final anal- 
ysis, due to its large deviation from the general trend of the 
data, which might be caused by its incomplete sampling.) 
To determine these times of minima, first we removed the 
effects of the intrinsic brightness variations from the light 
curves, and then fitted each outer transit and occultation 
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Figure 1. Example of a primary (upper panel) and a secondary 
{lower panel) shallow minimum to illustrate the states before 
(black triangles) and after (red crosses) the detrending of the 
minima. The dashed line is the fit used for the detrending (see 
Section [211. 



event individually with our newly developed simultaneous 
light curve solution code. Both the code and the complete 
light curve analysis are described in Sect. |4] 

The determined times of minima are listed in Tables [T] 
and [2] for the close and the wide pairs, respectively. 



3 ANALYSIS OF THE ECLIPSE TIMING 
VARIATIONS (ETV) 

3.1 The close binary 

In order to study the eclipse timing variations (ETV), the 
following linear ephemeris was calculated for the shallow 
minima: 



M/iVi_ 



[BJD] = 2 455 051.23625 -f- 0'?905677 x E, (1) 



where E is the cycle number. The corresponding ETV dia- 
gram is plotted in Fig. O 

We see a sinusoidal variation with a period identical 
to the eclipsing period of the wide system. There is also a 
smaller, long-term variation, that might either be part of a 
longer period variation, or represent a secular trend, as is 
the case with several close binary systems. First we analyse 
the periodic behaviour of the ETV, and then the possible 
secular (parabolic) term will also be discussed. 

3.1.1 Short-period variations: General remarks 

To detect further periodicities, a discrete Fourier transform 
was calculated for the ETV curve. The resulting amplitude 
spectrum shows that the odd harmonics of the fundamen- 
tal frequency are also present (see Fig. [3)l, while only the 
first even harmonic (i.e. 2/o) exists, and its amplitude is 
smaller than that of the 3/o and 5/o components. To check 
whether this structure is a consequence of the non-uniform 
sampling (i.e., the missing data during the deep eclipses, 
when the eclipse-events of the close pair cannot be observed, 
see Fig. |4] below), we calculated a simple circular light-time 
orbit solution (i.e. we first fitted a sine curve with the fun- 
damental frequency of the DFT spectrum). Sampling this 
solution at the locations (i.e. cycle numbers) of the observed 
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Table 1. Times of minima for the close pair 
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.0010 


II 
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II 
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II 
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I 
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I 
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I 
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.0010 


I 
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II 
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.0010 


II 


2455016. 


.8215 


0, 


.0010 


I 


2455121. 


.4274 


0, 


.0010 


II 


2455151, 


.3110 


0. 


.0010 


II 


2454985, 


.1226 


0. 


,0010 


I 


2455019. 


.5391 


0. 


.0010 


I 


2455121. 


.8792 


0, 


.0010 


I 


2455151, 


.7693 


0. 


,0010 


I 


2454985, 


.5771 


0. 


.0010 


II 


2455019. 


.9919 


0. 


.0010 


II 


2455122. 


.3325 


0, 


.0010 


II 


2455152, 


.2221 


0. 


,0010 


II 


2454986, 


.0275 


0. 


.0010 


I 


2455020. 


.4449 


0. 


.0010 


I 


2455122. 


.7838 


0, 


.0010 


I 


2455152, 


.6729 


0. 


,0010 


I 


2454986, 


.4819 


0. 


.0010 


II 


2455020. 


.8964 


0, 


.0010 


II 


2455123. 


.2389 


0, 


.0010 


II 


2455153, 


.1263 


0. 


,0010 


II 


2454986, 


.9334 


0. 


.0010 


I 


2455093. 


.3498 


0, 


.0010 


II 


2455124. 


.5962 


0, 


.0010 


I 


2455153, 


.5789 


0. 


,0010 


I 


2454987, 


.3857 


0. 


.0010 


II 


2455093. 


.8030 


0. 


.0010 


I 


2455125. 


.0489 


0, 


.0010 


II 


2455154, 


.0318 


0. 


.0010 


II 


2454987, 


.8395 


0. 


.0010 


I 


2455094. 


.2549 


0. 


.0010 


II 


2455125. 


.5019 


0, 


.0010 


I 


2455156, 


.7506 


0. 


,0010 


II 


2454988, 


.2928 


0. 


.0010 


II 


2455094. 


.7076 


0, 


.0010 


I 


2455125. 


.9538 


0, 


.0010 


II 


2455157, 


.2025 


0. 


.0010 


I 


2454988, 


.7449 


0. 


.0010 


I 


2455095. 


.1606 


0. 


.0010 


II 


2455126. 


.4076 


0, 


.0010 


I 


2455157, 


.6564 


0. 


,0010 


II 


2454989, 


.1967 


0. 


.0010 


II 


2455095. 


.6140 


0. 


.0010 


I 


2455126. 


.8605 


0, 


.0010 


II 


2455160, 


.3727 


0. 


.0010 


II 


2454989, 


.6504 


0. 


.0010 


I 


2455096. 


.0667 


0. 


.0010 


II 


2455127. 


.3131 


0, 


.0010 


I 


2455160, 


.8254 


0. 


,0010 


I 


2454990, 


.1049 


0. 


.0010 


II 


2455096. 


.5190 


0, 


.0010 


I 


2455127. 


.7654 


0, 


.0010 


II 


2455161, 


.2782 


0. 


.0010 


II 


2454990, 


.5561 


0. 


.0010 


I 


2455096. 


.9724 


0. 


.0010 


II 


2455128, 


.2180 


0, 


.0010 


I 


2455161, 


.7294 


0. 


,0010 


I 


2454991, 


.0091 


0. 


.0010 


II 


2455097. 


.4240 


0. 


.0010 


I 


2455128. 


.6723 


0, 


.0010 


II 


2455162, 


.1839 


0. 


.0010 


II 


2454991, 


.4619 


0. 


.0010 


I 


2455097. 


.8796 


0. 


.0010 


II 


2455129. 


.1244 


0, 


.0010 


I 


2455162, 


.6355 


0. 


,0010 


I 


2454991, 


.9159 


0. 


.0010 


II 


2455098. 


.3312 


0, 


.0010 


I 


2455129. 


.5755 


0, 


.0010 


II 


2455163, 


.0890 


0. 


.0010 


II 


2454992, 


.3675 


0. 


.0010 


I 


2455098. 


.7829 


0. 


.0010 


II 


2455130. 


.0306 


0, 


.0010 


I 


2455163, 


.5420 


0. 


,0010 


I 


2454992, 


.8209 


0. 


.0010 


II 


2455099. 


.2375 


0. 


.0010 


I 


2455130. 


.4822 


0, 


.0010 


II 


2455163, 


.9946 


0. 


.0010 


II 


2454993, 


.2726 


0. 


.0010 


I 


2455099. 


.6907 


0. 


.0010 


II 


2455130, 


.9330 


0, 


.0010 


I 


2455164, 


.4462 


0. 


.0010 


I 


2454993, 


.7237 


0. 


,0010 


II 


2455100. 


.1415 


0. 


.0010 


I 


2455131. 


.3871 


0, 


.0010 


II 


2455164, 


.8980 


0. 


.0010 


II 


2454994, 


.1785 


0, 


,0010 


I 


2455101. 


.0502 


0. 


.0010 


I 


2455131. 


.8408 


0, 


.0010 


I 


2455165, 


.3521 


0. 


,0010 


I 
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Table 1. Times of minima for the close pair (continued) 



BJD 








lype 


BJD 








±ype 


BJD 










BJD 








Tvnp 

±ype 


2455165, 


.8057 


0. 


.0010 


II 


2455196. 


.5973 


0. 


.0010 


II 


2455229. 


.2014 


0, 


.0010 


II 


2455262, 


.2589 


0. 


,0010 


I 


2455166, 


.2571 


0. 


.0010 


I 


2455197. 


.0506 


0. 


.0010 


I 


2455229. 


.6530 


0, 


.0010 


I 


2455262, 


.7143 


0. 


,0010 


II 


2455166, 


.7118 


0. 


.0010 


II 


2455197. 


.5016 


0. 


.0010 


II 


2455234. 


.1820 


0, 


.0010 


I 


2455263, 


.1643 


0. 


,0010 


I 


2455167, 


.1625 


0. 


.0010 


I 


2455197. 


.9563 


0, 


.0010 


I 


2455234. 


.6354 


0, 


.0010 


II 


2455263, 


.6172 


0. 


.0010 


II 


2455167, 


.6143 


0. 


.0010 


II 


2455198. 


.8642 


0, 


.0010 


I 


2455235. 


.0871 


0, 


.0010 


I 


2455264, 


.0693 


0. 


,0010 


I 


2455168, 


.0689 


0. 


.0010 


I 


2455199. 


.3154 


0. 


.0010 


II 


2455235. 


.5403 


0, 


.0010 


II 


2455264, 


.5227 


0. 


.0010 


II 


2455168, 


.5204 


0. 


.0010 


II 


2455199. 


.7682 


0. 


.0010 


I 


2455235. 


.9946 


0, 


.0010 


I 


2455264, 


.9751 


0. 


,0010 


I 


2455168, 


.9745 


0. 


.0010 


I 


2455200. 


.2232 


0, 


.0010 


II 


2455236. 


.4474 


0, 


.0010 


II 


2455265, 


.4305 


0. 


.0010 


II 


2455169, 


.4273 


0. 


.0010 


II 


2455200. 


.6731 


0, 


.0010 


I 


2455236. 


.8990 


0, 


.0010 


I 


2455265, 


.8805 


0. 


,0010 


I 


2455169, 


.8793 


0. 


.0010 


I 


2455201. 


.1280 


0. 


.0010 


II 


2455237. 


.3518 


0, 


.0010 


II 


2455266, 


.3335 


0. 


.0010 


II 


2455170, 


.3314 


0, 


.0010 


II 


2455201. 


.5798 


0. 


.0010 


I 


2455237. 


.8052 


0, 


.0010 


I 


2455266, 


.7867 


0. 


,0010 


I 


2455170, 


.7846 


0. 


.0010 


I 


2455202. 


.0321 


0. 


.0010 


II 


2455238. 


.7127 


0, 


.0010 


I 


2455267, 


.2403 


0. 


,0010 


II 


2455171, 


.2377 


0. 


.0010 


II 


2455202. 


.4852 


0. 


.0010 


I 


2455239. 


.1626 


0, 


.0010 


II 


2455267, 


.6921 


0. 


.0010 


I 


2455171, 


.6922 


0. 


.0010 


I 


2455202. 


.9387 


0. 


.0010 


II 


2455239. 


.6170 


0, 


.0010 


I 


2455268, 


.1446 


0. 


,0010 


II 


2455172, 


.1444 


0. 


.0010 


II 


2455205, 


.6558 


0, 


.0010 


II 


2455240. 


.0709 


0, 


.0010 


II 


2455268, 


.5972 


0. 


,0010 


I 


2455172, 


.5969 


0. 


.0010 


I 


2455206. 


.1085 


0, 


.0010 


I 


2455240. 


.5226 


0, 


.0010 


I 


2455269, 


.0530 


0. 


,0010 


II 


2455173, 


.0483 


0. 


.0010 


II 


2455206. 


.5624 


0, 


.0010 


II 


2455240. 


.9766 


0, 


.0010 


II 


2455269, 


.9571 


0. 


,0010 


II 


2455173, 


.5019 


0. 


.0010 


I 


2455207. 


.0141 


0. 


.0010 


I 


2455241. 


.4291 


0, 


.0010 


I 


2455270, 


.4094 


0. 


,0010 


I 


2455173, 


.9555 


0. 


.0010 


II 


2455207, 


.4666 


0. 


.0010 


II 


2455241. 


.8823 


0, 


.0010 


II 


2455270, 


.8620 


0. 


,0010 


II 


2455174, 


.4083 


0. 


.0010 


I 


2455207. 


.9207 


0, 


.0010 


I 


2455242. 


.3343 


0, 


.0010 


I 


2455271, 


.3125 


0. 


,0010 


I 


2455174, 


.8618 


0. 


.0010 


II 


2455208. 


.3747 


0. 


.0010 


II 


2455242. 


.7872 


0, 


.0010 


II 


2455274, 


.0311 


0. 


,0010 


I 


2455175, 


.3120 


0. 


.0010 


I 


2455208. 


.8256 


0. 


.0010 


I 


2455243. 


.2405 


0, 


.0010 


I 


2455274, 


.4849 


0. 


,0010 


II 


2455176, 


.2170 


0. 


.0010 


I 


2455209, 


.2773 


0. 


.0010 


II 


2455243. 


.6902 


0, 


.0010 


II 


2455274, 


.9361 


0. 


.0010 


I 


2455176, 


.6701 


0. 


.0010 


II 


2455209. 


.7301 


0. 


.0010 


I 


2455244. 


.1458 


0, 


.0010 


I 


2455277, 


.2007 


0. 


,0010 


II 


2455177, 


.1247 


0. 


.0010 


I 


2455210. 


.1846 


0, 


.0010 


II 


2455244. 


.5992 


0, 


.0010 


II 


2455278, 


.1082 


0. 


.0010 


II 


2455177, 


.5766 


0. 


.0010 


II 


2455210. 


.6363 


0, 


.0010 


I 


2455245, 


.0515 


0, 


.0010 


I 


2455279, 


.9181 


0. 


,0010 


II 


2455178, 


.0299 


0. 


.0010 


I 


2455211. 


.0904 


0. 


.0010 


II 


2455245. 


.5051 


0, 


.0010 


II 


2455280, 


.8252 


0. 


.0010 


II 


2455178, 


.4846 


0. 


.0010 


II 


2455211. 


.5409 


0. 


.0010 


I 


2455245. 


.9575 


0, 


.0010 


I 


2455281, 


.7304 


0. 


,0010 


II 


2455178, 


.9346 


0. 


.0010 


I 


2455211. 


.9941 


0. 


.0010 


II 


2455246. 


.4102 


0, 


.0010 


II 


2455282, 


.6374 


0. 


.0010 


II 


2455179, 


.3877 


0. 


.0010 


II 


2455212. 


.4471 


0. 


.0010 


I 


2455246. 


.8628 


0, 


.0010 


I 


2455283, 


.5419 


0. 


,0010 


II 


2455179, 


.8418 


0. 


.0010 


I 


2455212. 


.8994 


0. 


.0010 


II 


2455247. 


.3158 


0, 


.0010 


II 


2455284, 


.4490 


0. 


.0010 


II 


2455180, 


.2948 


0. 


.0010 


II 


2455213. 


.3535 


0. 


.0010 


I 


2455247. 


.7698 


0, 


.0010 


I 


2455285, 


.3552 


0. 


,0010 


II 


2455180, 


.7447 


0. 


.0010 


I 


2455213. 


.8046 


0, 


.0010 


II 


2455248. 


.2229 


0, 


.0010 


II 


2455286, 


.2593 


0. 


.0010 


II 


2455184, 


.8249 


0. 


.0010 


II 


2455214. 


.2586 


0, 


.0010 


I 


2455248. 


.6748 


0, 


.0010 


I 


2455287, 


.1656 


0. 


,0010 


II 


2455185, 


.2776 


0. 


.0010 


I 


2455214. 


.7105 


0. 


.0010 


II 


2455250. 


.9482 


0, 


.0010 


II 


2455288, 


.0726 


0. 


.0010 


II 


2455185, 


.7273 


0. 


.0010 


II 


2455215, 


.1632 


0. 


.0010 


I 


2455251. 


.3928 


0, 


.0010 


I 


2455288, 


.9782 


0. 


,0010 


II 


2455186, 


.1799 


0. 


.0010 


I 


2455215. 


,6139 


0, 


.0010 


II 


2455251. 


.8454 


0, 


.0010 


II 


2455289, 


.8833 


0. 


.0010 


II 


2455186, 


.6346 


0. 


.0010 


II 


2455216. 


.0699 


0, 


.0010 


I 


2455252. 


.2979 


0, 


.0010 


I 


2455290, 


.7894 


0. 


,0010 


II 


2455187, 


.0869 


0. 


.0010 


I 


2455217, 


.4272 


0. 


.0010 


II 


2455252. 


.7506 


0, 


.0010 


II 


2455291, 


.6945 


0. 


.0010 


II 


2455187, 


.5412 


0. 


.0010 


II 


2455217. 


.8792 


0. 


.0010 


I 


2455253. 


.2041 


0, 


.0010 


I 


2455292, 


.5995 


0. 


.0010 


II 


2455187, 


.9934 


0. 


.0010 


I 


2455218. 


.3336 


0, 


.0010 


II 


2455253. 


.6582 


0, 


.0010 


II 


2455293, 


.5069 


0. 


.0010 


II 


2455188, 


.4450 


0. 


.0010 


II 


2455218. 


.7867 


0, 


.0010 


I 


2455254. 


.1080 


0, 


.0010 


I 


2455297, 


.1281 


0. 


,0010 


II 


2455188, 


.8981 


0. 


.0010 


I 


2455219. 


.2412 


0. 


.0010 


II 


2455254. 


.5616 


0, 


.0010 


II 


2455298, 


.0349 


0. 


.0010 


II 


2455189, 


.3515 


0. 


.0010 


II 


2455219. 


.6907 


0. 


.0010 


I 


2455255. 


.0146 


0, 


.0010 


I 


2455298, 


.9417 


0. 


,0010 


II 


2455189, 


.8027 


0. 


.0010 


I 


2455220. 


.1447 


0, 


.0010 


II 


2455255. 


.4681 


0, 


.0010 


II 


2455299, 


.8474 


0. 


.0010 


II 


2455190, 


.2577 


0. 


.0010 


II 


2455220. 


.5967 


0, 


.0010 


I 


2455255. 


.9214 


0, 


.0010 


I 


2455300, 


.7532 


0. 


,0010 


II 


2455190, 


.7087 


0, 


.0010 


I 


2455221. 


.0499 


0. 


.0010 


II 


2455256. 


.3741 


0, 


.0010 


II 


2455301, 


.6557 


0. 


.0010 


II 


2455191, 


.1655 


0, 


.0010 


II 


2455221. 


.5019 


0. 


.0010 


I 


2455256. 


.8255 


0, 


.0010 


I 


2455302, 


.5648 


0. 


,0010 


II 


2455191, 


.6160 


0. 


,0010 


I 


2455221. 


.9583 


0, 


.0010 


II 


2455257. 


.2747 


0, 


.0010 


II 


2455303, 


.4691 


0. 


.0010 


II 


2455192, 


.0681 


0. 


.0010 


II 


2455222. 


.4097 


0, 


.0010 


I 


2455257. 


.7299 


0, 


.0010 


I 


2455305, 


.2800 


0. 


.0010 


II 


2455192, 


.5204 


0. 


.0010 


I 


2455222. 


.8609 


0. 


.0010 


II 


2455258. 


.1856 


0, 


.0010 


II 


2455306, 


.1858 


0. 


.0010 


II 


2455192, 


.9764 


0. 


.0010 


II 


2455223. 


.3136 


0. 


.0010 


I 


2455258. 


.6384 


0, 


.0010 


I 


2455307, 


.0881 


0. 


,0010 


II 


2455193, 


.4266 


0. 


.0010 


I 


2455223. 


.7642 


0, 


.0010 


II 


2455259. 


.0912 


0, 


.0010 


II 


2455309, 


.8044 


0. 


.0010 


II 


2455193, 


.8803 


0. 


.0010 


II 


2455224. 


.2188 


0. 


.0010 


I 


2455259. 


.5424 


0, 


.0010 


I 


2455310, 


.7133 


0. 


,0010 


II 


2455194, 


.3331 


0. 


.0010 


I 


2455224. 


.6721 


0. 


.0010 


II 


2455259. 


.9961 


0, 


.0010 


II 


2455311, 


.6198 


0. 


.0010 


II 


2455194, 


.7853 


0. 


.0010 


II 


2455225. 


.1257 


0. 


.0010 


I 


2455260. 


.4492 


0, 


.0010 


I 


2455313, 


.4275 


0. 


,0010 


II 


2455195, 


.2387 


0. 


.0010 


I 


2455226. 


.0314 


0, 


.0010 


I 


2455260. 


.8994 


0, 


.0010 


II 


2455314, 


.3352 


0. 


.0010 


II 


2455195, 


.6946 


0. 


.0010 


II 


2455228. 


.2937 


0, 


.0010 


II 


2455261. 


.3547 


0, 


.0010 


I 


2455315, 


.2384 


0. 


,0010 


II 


2455196, 


.1452 


0. 


.0010 


I 


2455228, 


.7472 


0. 


.0010 


I 


2455261. 


.8072 


0, 


.0010 


II 


2455316, 


.1422 


0. 


.0010 


II 
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Table 1. Times of minima for the close pair (continued) 



BJD 








lype 


BJD 








lype 


BJD 










BJD 








Tvnp 

lype 


2455317, 


.0475 


0. 


.0010 


II 


2455377. 


.7335 


0. 


.0010 


II 


2455408. 


.0720 


0, 


.0010 


I 


2455439, 


.7717 


0. 


,0010 


I 


2455319, 


.7697 


0. 


.0010 


II 


2455378, 


.1862 


0. 


.0010 


I 


2455408. 


.5241 


0, 


.0010 


II 


2455440, 


.2272 


0. 


.0010 


II 


2455320, 


.6756 


0. 


.0010 


II 


2455378. 


.6391 


0, 


.0010 


II 


2455411, 


.2440 


0, 


.0010 


II 


2455440, 


.6775 


0. 


,0010 


I 


2455321, 


.5787 


0. 


.0010 


II 


2455379. 


.0916 


0, 


.0010 


I 


2455411, 


.6947 


0, 


.0010 


I 


2455441, 


.1326 


0. 


,0010 


II 


2455322, 


.4873 


0. 


.0010 


II 


2455379. 


.5443 


0, 


.0010 


II 


2455412, 


.1468 


0, 


.0010 


II 


2455441, 


.5837 


0. 


,0010 


I 


2455323, 


.3916 


0. 


.0010 


II 


2455379. 


.9967 


0. 


.0010 


I 


2455412, 


.6000 


0, 


.0010 


I 


2455442, 


.0395 


0. 


,0010 


II 


2455324, 


.2975 


0. 


.0010 


II 


2455380. 


.4486 


0. 


.0010 


II 


2455413. 


.0556 


0, 


.0010 


II 


2455442, 


.4897 


0. 


,0010 


I 


2455325, 


.2031 


0. 


.0010 


II 


2455380. 


,9024 


0, 


.0010 


I 


2455413. 


.5075 


0, 


.0010 


I 


2455442, 


.9436 


0. 


,0010 


II 


2455326, 


.1094 


0. 


.0010 


II 


2455381. 


.3567 


0, 


.0010 


II 


2455413. 


.9629 


0, 


.0010 


II 


2455443, 


.3940 


0. 


.0010 


I 


2455327, 


.0152 


0. 


.0010 


II 


2455381. 


.8100 


0. 


.0010 


I 


2455414, 


.4120 


0, 


.0010 


I 


2455443, 


.8493 


0. 


.0010 


II 


2455327, 


.9208 


0, 


.0010 


II 


2455382, 


.2629 


0. 


.0010 


II 


2455414, 


.8646 


0, 


.0010 


II 


2455444, 


.3003 


0. 


.0010 


I 


2455328, 


.8275 


0. 


.0010 


II 


2455382, 


.7152 


0. 


.0010 


I 


2455415. 


.3159 


0, 


.0010 


I 


2455444, 


.7546 


0. 


,0010 


II 


2455329, 


.7327 


0. 


.0010 


II 


2455383. 


.1670 


0. 


.0010 


II 


2455415. 


.7715 


0, 


.0010 


II 


2455445, 


.2046 


0. 


,0010 


I 


2455330, 


.6358 


0. 


.0010 


II 


2455383. 


.6215 


0. 


.0010 


I 


2455416. 


.2234 


0, 


.0010 


I 


2455445, 


.6593 


0. 


,0010 


II 


2455331, 


.5422 


0. 


.0010 


II 


2455384, 


.0749 


0, 


.0010 


II 


2455416. 


.6769 


0, 


.0010 


II 


2455446, 


.1110 


0. 


,0010 


I 


2455332, 


.4492 


0. 


.0010 


II 


2455384, 


.5267 


0, 


.0010 


I 


2455417. 


.1288 


0, 


.0010 


I 


2455446, 


.5625 


0. 


,0010 


II 


2455333, 


.3539 


0. 


.0010 


II 


2455384, 


.9805 


0, 


.0010 


II 


2455417. 


.5823 


0, 


.0010 


II 


2455447, 


.0161 


0. 


,0010 


I 


2455334, 


.2614 


0. 


.0010 


II 


2455385. 


.4313 


0. 


.0010 


I 


2455418. 


.0353 


0, 


.0010 


I 


2455447, 


.4660 


0. 


,0010 


II 


2455335, 


.1677 


0. 


.0010 


II 


2455385. 


.8849 


0. 


.0010 


II 


2455418. 


.4868 


0, 


.0010 


II 


2455447, 


.9214 


0. 


,0010 


I 


2455336, 


.0729 


0. 


.0010 


II 


2455388. 


.6028 


0. 


.0010 


II 


2455418. 


.9404 


0, 


.0010 


I 


2455448, 


.3760 


0. 


,0010 


II 


2455337, 


.8848 


0. 


.0010 


II 


2455389. 


.0558 


0. 


.0010 


I 


2455419. 


.3938 


0, 


.0010 


II 


2455448, 


.8280 


0. 


.0010 


I 


2455338, 


.7889 


0. 


.0010 


II 


2455389. 


.5085 


0. 


.0010 


II 


2455419. 


.8471 


0, 


.0010 


I 


2455449, 


.2839 


0. 


,0010 


II 


2455339, 


.6974 


0. 


.0010 


II 


2455389. 


.9616 


0. 


.0010 


I 


2455420. 


.3012 


0, 


.0010 


II 


2455449, 


.7337 


0. 


.0010 


I 


2455342, 


.4122 


0. 


.0010 


II 


2455390. 


.4144 


0. 


.0010 


II 


2455420. 


.7525 


0, 


.0010 


I 


2455450, 


.1861 


0. 


,0010 


II 


2455343, 


.3198 


0. 


.0010 


II 


2455390. 


.8672 


0, 


.0010 


I 


2455421, 


.2057 


0, 


.0010 


II 


2455450, 


.6387 


0. 


.0010 


I 


2455344, 


.2237 


0. 


.0010 


II 


2455391. 


,3201 


0, 


.0010 


II 


2455421, 


.6579 


0, 


.0010 


I 


2455451, 


.0936 


0. 


,0010 


II 


2455345, 


.1306 


0. 


.0010 


II 


2455391. 


.7716 


0. 


.0010 


I 


2455422, 


.1103 


0, 


.0010 


II 


2455451, 


.5456 


0. 


.0010 


I 


2455346, 


.0355 


0. 


.0010 


II 


2455392, 


.2254 


0. 


.0010 


II 


2455422, 


.5636 


0, 


.0010 


I 


2455451, 


.9975 


0. 


,0010 


II 


2455346, 


.9414 


0. 


.0010 


II 


2455392, 


.6779 


0. 


.0010 


I 


2455423. 


.0166 


0, 


.0010 


II 


2455452, 


.4485 


0. 


.0010 


I 


2455347, 


.8465 


0. 


.0010 


II 


2455393. 


.1300 


0. 


.0010 


II 


2455423. 


.4695 


0, 


.0010 


I 


2455452, 


.9030 


0. 


.0010 


II 


2455348, 


.7514 


0. 


.0010 


II 


2455393. 


.5834 


0. 


.0010 


I 


2455423. 


.9224 


0, 


.0010 


II 


2455453, 


.3558 


0. 


.0010 


I 


2455349, 


.6589 


0. 


.0010 


II 


2455394, 


.0365 


0. 


.0010 


II 


2455424, 


.3746 


0, 


.0010 


I 


2455453, 


.8064 


0. 


,0010 


II 


2455350, 


.5622 


0. 


.0010 


II 


2455394, 


.4891 


0, 


.0010 


I 


2455424, 


.8276 


0, 


.0010 


II 


2455454, 


.2636 


0. 


.0010 


I 


2455351, 


.4699 


0. 


.0010 


II 


2455394, 


.9404 


0, 


.0010 


II 


2455425. 


.2819 


0, 


.0010 


I 


2455456, 


.5287 


0. 


,0010 


II 


2455352, 


.3738 


0. 


.0010 


II 


2455395. 


.3944 


0. 


.0010 


I 


2455425. 


.7354 


0, 


.0010 


II 


2455456, 


.9785 


0. 


.0010 


I 


2455353, 


.2807 


0. 


.0010 


II 


2455395. 


.8469 


0. 


.0010 


II 


2455426. 


.6409 


0, 


.0010 


II 


2455457, 


.4319 


0. 


,0010 


II 


2455354, 


.1831 


0. 


.0010 


II 


2455396. 


.2995 


0. 


.0010 


I 


2455427. 


.0935 


0, 


.0010 


I 


2455457, 


.8844 


0. 


,0010 


I 


2455355, 


.0909 


0. 


.0010 


II 


2455396. 


.7526 


0. 


.0010 


II 


2455427. 


.5467 


0, 


.0010 


II 


2455458, 


.3384 


0. 


,0010 


II 


2455355, 


.9974 


0. 


.0010 


II 


2455397, 


.2049 


0. 


.0010 


I 


2455427. 


.9997 


0, 


.0010 


I 


2455458, 


.7885 


0. 


,0010 


I 


2455356, 


.9013 


0. 


.0010 


II 


2455397. 


.6620 


0. 


.0010 


II 


2455428. 


.4520 


0, 


.0010 


II 


2455459, 


.2414 


0. 


,0010 


II 


2455357, 


.8074 


0. 


.0010 


II 


2455398. 


.1130 


0, 


.0010 


I 


2455428. 


.9047 


0, 


.0010 


I 


2455459, 


.6944 


0. 


,0010 


I 


2455358, 


.7129 


0. 


.0010 


II 


2455398. 


.5619 


0, 


.0010 


II 


2455429, 


.3569 


0, 


.0010 


II 


2455460, 


.1483 


0. 


,0010 


II 


2455359, 


.6172 


0. 


.0010 


II 


2455399. 


.0156 


0. 


.0010 


I 


2455429. 


.8087 


0, 


.0010 


I 


2455460, 


.6005 


0. 


,0010 


I 


2455360, 


.5266 


0. 


.0010 


II 


2455399. 


.4689 


0. 


.0010 


II 


2455430. 


.2654 


0, 


.0010 


II 


2455461, 


.0525 


0. 


,0010 


II 


2455361, 


.4287 


0. 


.0010 


II 


2455399. 


.9218 


0, 


.0010 


I 


2455430. 


.7160 


0, 


.0010 


I 


2455461, 


.5072 


0. 


,0010 


I 


2455362, 


.3342 


0. 


.0010 


II 


2455401. 


.7331 


0, 


.0010 


I 


2455431. 


.1699 


0, 


.0010 


II 


2455461, 


.9624 


0. 


.0010 


II 


2455365, 


.0530 


0, 


.0010 


II 


2455402, 


.1855 


0. 


.0010 


II 


2455433. 


.8849 


0, 


.0010 


II 


2455462, 


.4128 


0. 


,0010 


I 


2455365, 


.9592 


0, 


.0010 


II 


2455402, 


.6390 


0. 


.0010 


I 


2455434, 


.3388 


0, 


.0010 


I 


2455462, 


.8674 


0. 


,0010 


II 


2455366, 


.8637 


0. 


,0010 


II 


2455403. 


.0901 


0, 


.0010 


II 


2455434, 


.7910 


0, 


.0010 


II 


2455463, 


.7717 


0. 


,0005 


II 


2455371, 


.8469 


0. 


.0010 


I 


2455403. 


.5441 


0, 


.0010 


I 


2455435. 


.2442 


0, 


.0010 


I 


2455464, 


.2261 


0. 


,0005 


I 


2455372, 


.2998 


0. 


.0010 


II 


2455403. 


.9981 


0. 


.0010 


II 


2455435. 


.6987 


0, 


.0010 


II 


2455464, 


.6770 


0. 


,0005 


II 


2455374, 


.1114 


0. 


.0010 


II 


2455404. 


.4514 


0. 


.0010 


I 


2455436. 


.1498 


0, 


.0010 


I 


2455465, 


.1315 


0. 


,0005 


I 


2455374, 


.5610 


0. 


.0010 


I 


2455404, 


.9002 


0, 


.0010 


II 


2455436. 


.6029 


0, 


.0010 


II 


2455465, 


.5842 


0. 


.0005 


II 


2455375, 


.0148 


0. 


.0010 


II 


2455405. 


.3552 


0, 


.0010 


I 


2455437. 


.0560 


0, 


.0010 


I 


2455466, 


.0369 


0. 


,0005 


I 


2455375, 


.4687 


0. 


.0010 


I 


2455405, 


.8120 


0. 


.0010 


II 


2455437. 


.5097 


0, 


.0010 


II 


2455466, 


.4878 


0. 


.0005 


II 


2455375, 


.9223 


0. 


.0010 


II 


2455406, 


.2622 


0. 


.0010 


I 


2455437. 


.9610 


0, 


.0010 


I 


2455466, 


.9428 


0. 


,0005 


I 


2455376, 


.3742 


0. 


.0010 


I 


2455406. 


.7133 


0. 


.0010 


II 


2455438. 


.4152 


0, 


.0010 


II 


2455467, 


.3976 


0. 


.0005 


II 


2455376, 


.8286 


0. 


.0010 


II 


2455407. 


.1677 


0. 


.0010 


I 


2455438. 


.8682 


0, 


.0010 


I 


2455467, 


.8483 


0. 


,0005 


I 


2455377, 


.2789 


0. 


.0010 


I 


2455407. 


.6199 


0. 


.0010 


II 


2455439. 


.3209 


0, 


.0010 


II 


2455468, 


.3019 


0. 


.0005 


II 
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Table 1. Times of minima for the close pair (continued) 



BJD 








lype 


BJD 








±ype 


BJD 










BJD 








Tvnp 

±ype 


2455468, 


.7542 


0. 


.0005 


I 


2455501, 


.3575 


0. 


.0005 


I 


2455531. 


.2464 


0. 


,0005 


I 


2455578, 


.3420 


0. 


,0005 


I 


2455469, 


.2073 


0. 


.0005 


II 


2455501. 


.8090 


0. 


.0005 


II 


2455531. 


.6999 


0. 


,0005 


II 


2455578, 


.7932 


0. 


,0005 


II 


2455469, 


.6603 


0. 


.0005 


I 


2455502. 


.2632 


0. 


.0005 


I 


2455532. 


.1515 


0. 


,0005 


I 


2455579, 


.6986 


0. 


.0005 


II 


2455470, 


.1139 


0. 


.0005 


II 


2455502. 


.7161 


0, 


.0005 


II 


2455532. 


.6062 


0. 


,0005 


II 


2455580, 


.1521 


0. 


,0005 


I 


2455470, 


.5665 


0. 


.0005 


I 


2455503. 


.1681 


0, 


.0005 


I 


2455533. 


.0574 


0. 


,0005 


I 


2455580, 


.6049 


0. 


,0005 


II 


2455471, 


.0186 


0. 


.0005 


II 


2455503. 


.6213 


0. 


.0005 


II 


2455533. 


.5100 


0. 


,0005 


II 


2455581, 


.0574 


0. 


,0005 


I 


2455471, 


.4718 


0. 


.0005 


I 


2455504. 


.0741 


0. 


.0005 


I 


2455533. 


.9625 


0. 


,0005 


I 


2455581, 


.5105 


0. 


,0005 


II 


2455471, 


.9251 


0. 


,0005 


II 


2455504. 


.5265 


0, 


.0005 


II 


2455534. 


.4158 


0. 


,0005 


II 


2455581, 


.9643 


0. 


,0005 


I 


2455472, 


.3781 


0. 


.0005 


I 


2455504. 


.9801 


0, 


.0005 


I 


2455534. 


.8681 


0. 


,0005 


I 


2455582, 


.4163 


0. 


,0005 


II 


2455472, 


.8314 


0. 


.0005 


II 


2455505. 


.4323 


0. 


.0005 


II 


2455535. 


.3218 


0. 


,0005 


II 


2455582, 


.8691 


0. 


,0005 


I 


2455473, 


.2831 


0, 


.0005 


I 


2455505. 


.8865 


0. 


.0005 


I 


2455535. 


.7741 


0. 


,0005 


I 


2455583, 


.3218 


0. 


,0005 


II 


2455473, 


.7377 


0. 


,0005 


II 


2455506. 


.3385 


0. 


.0005 


II 


2455536. 


.2268 


0. 


,0005 


II 


2455583, 


.7738 


0. 


,0005 


I 


2455474, 


.6430 


0. 


.0005 


II 


2455506. 


.7915 


0. 


.0005 


I 


2455536. 


.6788 


0. 


,0005 


I 


2455584, 


.2276 


0. 


,0005 


II 


2455475, 


.0950 


0. 


.0005 


I 


2455507. 


.2443 


0. 


.0005 


II 


2455537. 


.1329 


0. 


,0005 


II 


2455584, 


.6789 


0. 


.0005 


I 


2455475, 


.5475 


0. 


.0005 


II 


2455507. 


.6974 


0, 


.0005 


I 


2455537. 


.5842 


0. 


,0005 


I 


2455585, 


.1325 


0. 


,0005 


II 


2455476, 


.0015 


0. 


.0005 


I 


2455508. 


.6032 


0, 


.0005 


I 


2455538. 


.0368 


0. 


,0005 


II 


2455585, 


.5851 


0. 


.0005 


I 


2455478, 


.7173 


0. 


.0005 


I 


2455509. 


.0573 


0, 


.0005 


II 


2455538. 


.4896 


0. 


,0005 


I 


2455586, 


.4906 


0. 


,0005 


I 


2455479, 


.1718 


0. 


.0005 


II 


2455509. 


.5096 


0. 


.0005 


I 


2455538. 


.9423 


0. 


,0005 


II 


2455586, 


.9434 


0. 


.0005 


II 


2455479, 


.6237 


0, 


.0005 


I 


2455509. 


.9625 


0. 


.0005 


II 


2455539. 


.3957 


0. 


,0005 


I 


2455587, 


.3954 


0. 


,0005 


I 


2455480, 


.0770 


0. 


.0005 


II 


2455510. 


.4144 


0. 


.0005 


I 


2455539. 


.8489 


0. 


,0005 


II 


2455587, 


.8488 


0. 


.0005 


II 


2455480, 


.5302 


0. 


,0005 


I 


2455510. 


.8680 


0. 


.0005 


II 


2455540. 


.3005 


0. 


,0005 


I 


2455588, 


.3015 


0. 


,0005 


I 


2455480, 


.9821 


0. 


.0005 


II 


2455511. 


.3213 


0, 


.0005 


I 


2455540. 


.7539 


0. 


,0005 


II 


2455588, 


.7542 


0. 


.0005 


II 


2455481, 


.4349 


0. 


.0005 


I 


2455511. 


.7738 


0. 


.0005 


II 


2455541. 


.2068 


0. 


,0005 


I 


2455589, 


.2072 


0. 


,0005 


I 


2455481, 


.8883 


0. 


.0005 


II 


2455512, 


.2263 


0. 


.0005 


I 


2455541. 


.6590 


0. 


,0005 


II 


2455589, 


.6597 


0. 


.0005 


II 


2455482, 


.3405 


0. 


.0005 


I 


2455512. 


,6809 


0, 


.0005 


II 


2455542. 


.1120 


0. 


,0005 


I 


2455592, 


.3773 


0. 


,0005 


II 


2455482, 


.7922 


0. 


.0005 


II 


2455513. 


.1336 


0, 


.0005 


I 


2455542. 


.5656 


0. 


,0005 


II 


2455592, 


.8297 


0. 


.0005 


I 


2455483, 


.2465 


0. 


.0005 


I 


2455513. 


.5859 


0. 


.0005 


II 


2455543. 


.0184 


0. 


,0005 


I 


2455593, 


.2826 


0. 


,0005 


II 


2455483, 


.7020 


0. 


.0005 


II 


2455514. 


.0381 


0. 


.0005 


I 


2455543. 


.4702 


0. 


,0005 


II 


2455593, 


.7360 


0. 


.0005 


I 


2455484, 


.1517 


0. 


.0005 


I 


2455514. 


.4902 


0, 


.0005 


II 


2455547. 


.0937 


0. 


,0005 


II 


2455596, 


.9066 


0. 


.0005 


II 


2455484, 


.6053 


0. 


.0005 


II 


2455514. 


.9437 


0, 


.0005 


I 


2455547. 


.5458 


0. 


,0005 


I 


2455597, 


.3587 


0. 


,0005 


I 


2455485, 


.0578 


0. 


.0005 


I 


2455515. 


.3971 


0. 


.0005 


II 


2455547. 


.9990 


0. 


,0005 


II 


2455597, 


.8120 


0. 


.0005 


II 


2455485, 


.5097 


0. 


.0005 


II 


2455515. 


.8504 


0. 


.0005 


I 


2455548. 


.4525 


0. 


,0005 


I 


2455598, 


.2641 


0. 


,0005 


I 


2455485, 


.9627 


0. 


.0005 


I 


2455516. 


.3035 


0, 


.0005 


II 


2455548. 


.9056 


0. 


,0005 


II 


2455598, 


.7184 


0. 


.0005 


II 


2455486, 


.4147 


0. 


.0005 


II 


2455516. 


.7557 


0, 


.0005 


I 


2455549. 


.3570 


0. 


,0005 


I 


2455599, 


.1709 


0. 


,0005 


I 


2455486, 


.8686 


0. 


.0005 


I 


2455517, 


.2081 


0. 


.0005 


II 


2455549. 


.8100 


0. 


,0005 


II 


2455599, 


.6232 


0. 


.0005 


II 


2455487, 


.3220 


0. 


.0005 


II 


2455517. 


.6619 


0. 


.0005 


I 


2455550. 


.2637 


0. 


,0005 


I 


2455600, 


.0759 


0. 


,0005 


I 


2455487, 


.7737 


0. 


,0005 


I 


2455518. 


.1156 


0. 


.0005 


II 


2455550. 


.7169 


0. 


,0005 


II 


2455600, 


.5298 


0. 


.0005 


II 


2455488, 


.2261 


0. 


.0005 


II 


2455518. 


.5680 


0. 


.0005 


I 


2455551. 


.1691 


0. 


,0005 


I 


2455600, 


.9820 


0. 


.0005 


I 


2455488, 


.6789 


0. 


.0005 


I 


2455519. 


.0206 


0. 


.0005 


II 


2455551. 


.6238 


0. 


,0005 


II 


2455601, 


.4349 


0. 


.0005 


II 


2455489, 


.1323 


0. 


.0005 


II 


2455519. 


.4744 


0. 


.0005 


I 


2455552. 


.0754 


0. 


,0005 


I 


2455601, 


.8876 


0. 


,0005 


I 


2455489, 


.5846 


0. 


.0005 


I 


2455519. 


.9260 


0. 


.0005 


II 


2455569. 


.7379 


0. 


,0005 


II 


2455602, 


.3424 


0. 


.0005 


II 


2455490, 


.0371 


0. 


.0005 


II 


2455520. 


.3791 


0. 


.0005 


I 


2455570. 


.1922 


0. 


,0005 


I 


2455602, 


.7933 


0. 


,0005 


I 


2455490, 


.4899 


0. 


.0005 


I 


2455520. 


.8313 


0. 


.0005 


II 


2455570. 


.6430 


0. 


,0005 


II 


2455603, 


.2462 


0. 


.0005 


II 


2455490, 


.9432 


0. 


.0005 


II 


2455521. 


.2845 


0. 


.0005 


I 


2455571. 


.0980 


0. 


,0005 


I 


2455603, 


.6989 


0. 


,0005 


I 


2455491, 


.3953 


0. 


.0005 


I 


2455524. 


.4545 


0. 


.0005 


II 


2455571. 


.5494 


0. 


,0005 


II 


2455604, 


.1543 


0. 


.0005 


II 


2455491, 


.8476 


0. 


.0005 


II 


2455524. 


.9077 


0. 


.0005 


I 


2455572. 


.0033 


0. 


,0005 


I 


2455604, 


.6052 


0. 


.0005 


I 


2455492, 


.3012 


0, 


.0005 


I 


2455525. 


.3594 


0. 


.0005 


II 


2455572. 


.4559 


0. 


,0005 


II 


2455605, 


.0599 


0. 


.0005 


II 


2455492, 


.7532 


0, 


.0005 


II 


2455525. 


.8125 


0. 


.0005 


I 


2455572. 


.9093 


0. 


,0005 


I 


2455605, 


.5115 


0. 


,0005 


I 


2455493, 


.2069 


0. 


.0005 


I 


2455526. 


.2669 


0, 


.0005 


II 


2455573. 


.3596 


0. 


,0005 


II 


2455605, 


.9647 


0. 


.0005 


II 


2455494, 


.5648 


0. 


.0005 


II 


2455526. 


.7179 


0, 


.0005 


I 


2455573. 


.8147 


0. 


,0005 


I 


2455606, 


.4173 


0. 


,0005 


I 


2455495, 


.0172 


0. 


.0005 


I 


2455527. 


.1722 


0. 


.0005 


II 


2455574. 


.2668 


0. 


,0005 


II 


2455606, 


.8715 


0. 


.0005 


II 


2455495, 


.4689 


0. 


.0005 


II 


2455527. 


.6246 


0. 


.0005 


I 


2455574. 


.7207 


0. 


,0005 


I 


2455607, 


.3236 


0. 


,0005 


I 


2455495, 


.9226 


0. 


.0005 


I 


2455528. 


.0772 


0. 


.0005 


II 


2455575. 


.1715 


0. 


,0005 


II 


2455607, 


.7765 


0. 


.0005 


II 


2455496, 


.3746 


0. 


.0005 


II 


2455528. 


.5289 


0. 


.0005 


I 


2455575. 


.6251 


0. 


,0005 


I 


2455608, 


.2291 


0. 


,0005 


I 


2455496, 


.8289 


0. 


.0005 


I 


2455528, 


.9816 


0. 


.0005 


II 


2455576. 


.0767 


0. 


,0005 


II 


2455608, 


.6823 


0. 


.0005 


II 


2455497, 


.2808 


0. 


.0005 


II 


2455529. 


.4345 


0. 


.0005 


I 


2455576. 


.5308 


0. 


,0005 


I 


2455609, 


.1350 


0. 


.0005 


I 


2455497, 


.7340 


0. 


.0005 


I 


2455529. 


.8870 


0, 


.0005 


II 


2455576. 


.9824 


0. 


,0005 


II 


2455609, 


.5880 


0. 


.0005 


II 


2455498, 


.1876 


0. 


.0005 


II 


2455530. 


.3409 


0. 


.0005 


I 


2455577. 


.4365 


0. 


,0005 


I 


2455610, 


.0405 


0. 


,0005 


I 


2455498, 


.6395 


0, 


.0005 


I 


2455530, 


.7942 


0. 


.0005 


II 


2455577. 


.8889 


0. 


,0005 


II 


2455610, 


.4941 


0. 


.0005 


II 
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Table 1. Times of minima for the close pair (continued) 



BJD 








lype 


BJD 








lype 


BJD 










BJD 








Tvnp 

lype 


2455610, 


.9466 


0. 


.0005 


I 


2455647, 


.1719 


0. 


.0005 


I 


2455677. 


.0593 


0. 


,0005 


I 


2455710, 


.5712 


0. 


,0005 


I 


2455611, 


.3995 


0. 


.0005 


II 


2455647, 


.6256 


0. 


.0005 


II 


2455677. 


.5115 


0. 


,0005 


II 


2455711, 


.0220 


0. 


,0005 


II 


2455611, 


.8529 


0. 


.0005 


I 


2455648. 


.0790 


0. 


.0005 


I 


2455678. 


.8702 


0. 


,0005 


I 


2455711, 


.4773 


0. 


,0005 


I 


2455612, 


.3055 


0. 


.0005 


II 


2455648. 


.5305 


0, 


.0005 


II 


2455679. 


.3227 


0. 


,0005 


II 


2455711, 


.9281 


0. 


,0005 


II 


2455615, 


.0236 


0. 


.0005 


II 


2455648. 


.9834 


0, 


.0005 


I 


2455679. 


.7759 


0. 


,0005 


I 


2455712, 


.3825 


0. 


,0005 


I 


2455615, 


.4750 


0. 


.0005 


I 


2455649. 


.4363 


0. 


.0005 


II 


2455680. 


.2286 


0. 


,0005 


II 


2455712, 


.8348 


0. 


.0005 


II 


2455615, 


.9271 


0. 


.0005 


II 


2455649. 


.8897 


0. 


.0005 


I 


2455680. 


.6813 


0. 


,0005 


I 


2455713, 


.2877 


0. 


,0005 


I 


2455616, 


.3809 


0. 


,0005 


I 


2455650. 


.3425 


0, 


.0005 


II 


2455683. 


.3990 


0. 


,0005 


I 


2455713, 


.7399 


0. 


,0005 


II 


2455617, 


.2863 


0. 


.0005 


I 


2455650. 


.7951 


0, 


.0005 


I 


2455683. 


.8511 


0. 


,0005 


II 


2455714, 


.1931 


0. 


.0005 


I 


2455617, 


.7385 


0. 


.0005 


II 


2455651. 


.2478 


0. 


.0005 


II 


2455684. 


.3045 


0. 


,0005 


I 


2455714, 


.6455 


0. 


,0005 


II 


2455618, 


.1924 


0, 


.0005 


I 


2455651. 


.7015 


0. 


.0005 


I 


2455684. 


.7566 


0. 


,0005 


II 


2455715, 


.0986 


0. 


.0005 


I 


2455618, 


.6455 


0. 


.0005 


II 


2455652. 


.1542 


0. 


.0005 


II 


2455685. 


.2096 


0. 


,0005 


I 


2455715, 


.5518 


0. 


,0005 


II 


2455619, 


.0977 


0. 


.0005 


I 


2455653, 


.0602 


0. 


.0005 


II 


2455685. 


.6613 


0. 


,0005 


II 


2455716, 


.4582 


0. 


.0005 


II 


2455619, 


.5523 


0. 


.0005 


II 


2455653, 


.5132 


0. 


.0005 


I 


2455686. 


.1165 


0. 


,0005 


I 


2455716, 


.9101 


0. 


,0005 


I 


2455620, 


.0035 


0. 


.0005 


I 


2455653. 


.9665 


0, 


.0005 


II 


2455686. 


.5683 


0. 


,0005 


II 


2455717, 


.3626 


0. 


,0005 


II 


2455620, 


.4562 


0. 


.0005 


II 


2455654. 


.4189 


0, 


.0005 


I 


2455687. 


.0217 


0. 


,0005 


I 


2455717, 


.8152 


0. 


,0005 


I 


2455620, 


.9085 


0. 


.0005 


I 


2455654. 


.8719 


0, 


.0005 


II 


2455687. 


.4731 


0. 


,0005 


II 


2455718, 


.2668 


0. 


.0005 


II 


2455621, 


.3623 


0. 


.0005 


II 


2455655. 


.3245 


0. 


.0005 


I 


2455687. 


.9278 


0. 


,0005 


I 


2455718, 


.7203 


0. 


,0005 


I 


2455621, 


.8147 


0, 


.0005 


I 


2455655. 


.7771 


0. 


.0005 


II 


2455688. 


.3801 


0. 


,0005 


II 


2455719, 


.1725 


0. 


.0005 


II 


2455622, 


.2670 


0. 


.0005 


II 


2455656. 


.2302 


0. 


.0005 


I 


2455688. 


.8330 


0. 


,0005 


I 


2455719, 


.6261 


0. 


,0005 


I 


2455622, 


.7196 


0. 


,0005 


I 


2455656. 


.6842 


0. 


.0005 


II 


2455689. 


.2853 


0. 


,0005 


II 


2455720, 


.0793 


0. 


.0005 


II 


2455623, 


.1728 


0. 


.0005 


II 


2455657. 


.1370 


0. 


.0005 


I 


2455689. 


.7390 


0. 


,0005 


I 


2455720, 


.5319 


0. 


,0005 


I 


2455623, 


.6254 


0. 


.0005 


I 


2455657, 


.5885 


0. 


.0005 


II 


2455690. 


.1923 


0. 


,0005 


II 


2455720, 


.9837 


0. 


.0005 


II 


2455624, 


.0789 


0. 


.0005 


II 


2455660. 


.7584 


0. 


.0005 


I 


2455690. 


.6444 


0. 


,0005 


I 


2455721, 


.4378 


0. 


,0005 


I 


2455624, 


.5310 


0. 


.0005 


I 


2455661. 


.2120 


0. 


.0005 


II 


2455691. 


.0981 


0. 


,0005 


II 


2455721, 


.8891 


0. 


.0005 


II 


2455624, 


.9824 


0. 


.0005 


II 


2455661. 


.6649 


0, 


.0005 


I 


2455691, 


.5507 


0. 


,0005 


I 


2455722, 


.3431 


0. 


,0005 


I 


2455625, 


.4362 


0. 


.0005 


I 


2455662. 


.1179 


0. 


.0005 


II 


2455692. 


.0032 


0. 


,0005 


II 


2455722, 


.7959 


0. 


.0005 


II 


2455625, 


.8874 


0. 


.0005 


II 


2455662, 


.5705 


0. 


.0005 


I 


2455692. 


.4558 


0. 


,0005 


I 


2455723, 


.2483 


0. 


,0005 


I 


2455626, 


.3419 


0. 


.0005 


I 


2455663. 


.0229 


0, 


.0005 


II 


2455692. 


.9090 


0. 


,0005 


II 


2455723, 


.7019 


0. 


.0005 


II 


2455626, 


.7942 


0. 


.0005 


II 


2455663. 


.4762 


0, 


.0005 


I 


2455693. 


.3622 


0. 


,0005 


I 


2455724, 


.1541 


0. 


,0005 


I 


2455627, 


.2478 


0. 


.0005 


I 


2455663. 


.9279 


0. 


.0005 


II 


2455693. 


.8154 


0. 


,0005 


II 


2455724, 


.6072 


0. 


.0005 


II 


2455627, 


.6989 


0. 


.0005 


II 


2455664. 


.3806 


0. 


.0005 


I 


2455694. 


.2680 


0. 


,0005 


I 


2455725, 


.0599 


0. 


,0005 


I 


2455628, 


.1529 


0. 


.0005 


I 


2455664. 


.8344 


0. 


.0005 


II 


2455694. 


.7206 


0. 


,0005 


II 


2455725, 


.5131 


0. 


.0005 


II 


2455628, 


.6048 


0. 


.0005 


II 


2455665. 


.2877 


0. 


.0005 


I 


2455695. 


.1733 


0. 


,0005 


I 


2455725, 


.9651 


0. 


,0005 


I 


2455629, 


.0582 


0. 


.0005 


I 


2455665. 


.7389 


0. 


.0005 


II 


2455695. 


.6255 


0. 


,0005 


II 


2455728, 


.6829 


0. 


.0005 


I 


2455629, 


.5099 


0. 


.0005 


II 


2455666. 


.1933 


0. 


.0005 


I 


2455696. 


.0795 


0. 


,0005 


I 


2455729, 


.1358 


0. 


,0005 


II 


2455629, 


.9642 


0. 


.0005 


I 


2455666. 


.6453 


0. 


.0005 


II 


2455696. 


.5309 


0. 


,0005 


II 


2455729, 


.5884 


0. 


.0005 


I 


2455630, 


.4169 


0. 


.0005 


II 


2455667. 


.0988 


0. 


.0005 


I 


2455696. 


.9853 


0. 


,0005 


I 


2455730, 


.0397 


0. 


,0005 


II 


2455630, 


.8700 


0. 


.0005 


I 


2455667, 


.5509 


0. 


.0005 


II 


2455697. 


.4391 


0. 


,0005 


II 


2455730, 


.4938 


0. 


.0005 


I 


2455631, 


.3218 


0. 


.0005 


II 


2455668. 


.0040 


0. 


.0005 


I 


2455697. 


.8911 


0. 


,0005 


I 


2455730, 


.9468 


0. 


,0005 


II 


2455631, 


.7755 


0. 


.0005 


I 


2455668. 


.4563 


0. 


.0005 


II 


2455698. 


.3443 


0. 


,0005 


II 


2455731, 


.3993 


0. 


.0005 


I 


2455632, 


.2287 


0. 


.0005 


II 


2455668. 


.9087 


0. 


.0005 


I 


2455698. 


.7968 


0. 


,0005 


I 


2455731, 


.8516 


0. 


,0005 


II 


2455632, 


.6800 


0. 


.0005 


I 


2455669. 


.3622 


0. 


.0005 


II 


2455699. 


.2494 


0. 


,0005 


II 


2455732, 


.3060 


0. 


.0005 


I 


2455633, 


.1334 


0. 


.0005 


II 


2455669. 


.8148 


0. 


.0005 


I 


2455699. 


.7030 


0. 


,0005 


I 


2455732, 


.7564 


0. 


,0005 


II 


2455633, 


.5857 


0. 


.0005 


I 


2455670. 


.2678 


0, 


.0005 


II 


2455700. 


.1542 


0. 


,0005 


II 


2455733, 


.2114 


0. 


.0005 


I 


2455634, 


.0386 


0. 


.0005 


II 


2455670. 


.7202 


0, 


.0005 


I 


2455700. 


.6085 


0. 


,0005 


I 


2455733, 


.6647 


0. 


,0005 


II 


2455634, 


.4916 


0, 


.0005 


I 


2455671. 


.1723 


0. 


.0005 


II 


2455701. 


.0623 


0. 


,0005 


II 


2455734, 


.1170 


0. 


.0005 


I 


2455641, 


.7369 


0, 


.0005 


I 


2455671, 


.6268 


0. 


.0005 


I 


2455701. 


.5145 


0. 


,0005 


I 


2455734, 


.5686 


0. 


,0005 


II 


2455642, 


.1894 


0. 


,0005 


II 


2455672. 


.0785 


0. 


.0005 


II 


2455701. 


.9651 


0. 


,0005 


II 


2455735, 


.0218 


0. 


,0005 


I 


2455642, 


.6431 


0. 


.0005 


I 


2455672. 


.5312 


0. 


.0005 


I 


2455702. 


.4201 


0. 


,0005 


I 


2455735, 


.4758 


0. 


,0005 


II 


2455643, 


.0961 


0. 


.0005 


II 


2455672. 


.9824 


0. 


.0005 


II 


2455702. 


.8724 


0. 


,0005 


II 


2455735, 


.9284 


0. 


,0005 


I 


2455643, 


.5483 


0. 


.0005 


I 


2455673. 


.4376 


0. 


.0005 


I 


2455703. 


.3256 


0. 


,0005 


I 


2455736, 


.3806 


0. 


,0005 


II 


2455644, 


.0018 


0. 


.0005 


II 


2455673. 


.8893 


0, 


.0005 


II 


2455706. 


.0422 


0. 


,0005 


I 


2455737, 


.2865 


0. 


,0005 


II 


2455644, 


.4544 


0. 


.0005 


I 


2455674. 


.3431 


0. 


.0005 


I 


2455706. 


.4951 


0. 


,0005 


II 


2455737, 


.7402 


0. 


.0005 


I 


2455644, 


.9086 


0. 


.0005 


II 


2455674. 


.7954 


0. 


.0005 


II 


2455706. 


.9482 


0. 


,0005 


I 


2455738, 


.1930 


0. 


.0005 


II 


2455645, 


.3605 


0. 


.0005 


I 


2455675. 


.2484 


0. 


.0005 


I 


2455708. 


.3071 


0. 


,0005 


II 


2455738, 


.6457 


0. 


.0005 


I 


2455645, 


.8139 


0. 


.0005 


II 


2455675. 


.6992 


0. 


.0005 


II 


2455708. 


.7601 


0. 


,0005 


I 












2455646, 


.2670 


0. 


.0005 


I 


2455676. 


.1537 


0. 


.0005 


I 


2455709. 


.2132 


0. 


,0005 


II 












2455646, 


.7209 


0, 


.0005 


II 


2455676. 


.6071 


0. 


.0005 


II 


2455709. 


.6667 


0. 


,0005 


I 













© 2012 RAS, MNRAS 000.[m23l 



Dynamical masses, absolute radii and 3D orbits of HD 181068 9 



0.006 



0.004 



Q1-6 (LC) 
Q7-9 (SC) 
Q1-9 fit 
secular trend 



0.002 - 




-0.002 



-0.004 



MIN| = 2455051 .23625 + 0.905677 E 



-100 



100 



200 



300 400 
Cycle Number 



500 



600 



700 



800 



Figure 2. Eclipse timing variations in the shallow minima. Triangles and circles mark the LC and SC data, respectively. The solid line 
stands for the Ql — Q9 fit, while the dashed line represents the weak secular (parabolic) trend. 



Table 2. Times of minima for the wide system 



BJD 


Cycle number" 


BJD 


Cycle number" 


2454977.0831 


-11.5 


2455363.5693 


-3.0 


2455022.5375 


-10.5 


2455386.3163 


-2.5 


2455045.2970 


-10.0 


2455409.0662 


-2.0 


2455068.0335 


-9.5 


2455431.7818 


-1.5 


2455113.5169 


-8.5 


2455454.5345 


-1.0 


2455136.2170 


-8.0 


2455477.2681 


-0.5 


2455158.9550 


-7.5 


2455499.9950 


0.0 


2455204.4405 


-6.5 


2455545.4559 


1.0 


2455227.1669 


-6.0 


2455590.9390 


2.0 


2455249.9048 


-5.5 


2455613.6734 


2.5 


2455272.6355 


-5.0 


2455659.1425 


3.5 


2455295.3893 


-4.5 


2455681.8955 


4.0 


2455318.1113 


-4.0 


2455704.6063 


4.5 


2455340.8384 


-3.5 


2455727.3559 


5.0 



° : half-integer values refer to secondary minima 

data, and calculating the DFT spectrum of this dataset, we 
found that the two spectra have very similar structure (see 
Fig. [3l , confirming our conjecture that the odd peaks are a 
data-sampling effect. Consequently, we restrict our analysis 
on the main peak (/o) and its second harmonic (2/o). 

Considering the fundamental term, it is clear that its 
main source should be the gravitational interaction between 
the inner, close binary, and the wider, more massive gi- 
ant star. This interaction has at least two consequences: 
(i) the geometrical light-time effect (LfTE), and (ii) a dy- 
namical effect, due to the gravitational perturbations of the 
third body on the close, inner binary. In the case of LITE, 
the amplitude of the effect increases wit h the separation , 
as seen in dozens of systems (see e. g. iQian et al.l |2012| : 
IPop fc Vamo3l2012l . for most recent examples). Conversely, 
the amplitudes of the dynamical terms scale with (P1/P2) 
which, due to various observational biases, makes this phe- 
nomenon difficult to detect with traditional ground-based 
o bservations. A de tailed analysis of this topic can be found 
in lBorkovits et all (,2003. , .2011 ). To our knowledge, the only 
system in which the dynamical effect was clearly detected by 
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Figure 3. The DFT amplitude spectrum of the ETV curve (lower 
solid line). In order to illustrate the possible data-sampling origin 
of the odd harmonics the spectrum of a similarly sampled sine 
function with /o-frequency is also plotted (upper dashed line). 



classical ground-based, small-aperturo photometric observa- 
tions, is lU Aurigae (Mayor 1990; Ozdcmir ct al. 2003). Nev- 
ertheless, for compact syste ms like the recently discovere d 
KOI- 126 (jCarter et al.ll201ll ). KOI-928 (jSteffen et al.ll201lh . 
the amplitude r atio may be reversed , as it was clearly shown 
for KOI-928 bylstefeneTaO (|201ll ). 

For HD 181068, we first consider the LITE contribution. 
Its shape and amplitude are: 



^JTVlite 

^LITE 



aa smz2 



ei) sinwB 



C 1-1-62 cos V2 

Irl X 10 — jj^smi2P2 (1 — e2Cos ojej 



2/3 



(2) 
(3) 



where ae, 12, 62, tJe, P2 are the semi-major axis, inclination, 
eccentricity, argument of periastron, and period of the bi- 
nary's orbit around the common centre of mass of the triple 
system. Furthermore, V2 is the true anomaly of the eclipsing 
pair in this orbit, ub = V2 + lob is its true longitude mea- 
sured from the intersection of the orbital plane and the plane 
of the sky, and c is the speed of light. (Inclination, eccen- 
tricity, period and true anomaly are simply given subscript 
2, because their values are identical to those of the rela- 
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tive wider orbit, traditionally centered on the inner binary.) 
Note also that in Eq. (|3} masses should be given in solar 
m asses, while per i od in days. Substituting the values found 
by iDerekas et al.l (|201ll ) (i.e., mx ~ 3Mq, tuab ~ 4.6M0, 



«2 



; 87?6, P2 « 45'?5, 62 = 0), we . 



Alite ~ 1.5 X 10" 



(4) 



or ~ 2.4 minutes. 

Now, considering the dynamical perturbation term, 
whose amplitude should be proportional to 



A 



1 ruA (1 - el 



,1/2 



27r mAB P2 (1 - elf'^ ' 
jBorkovits et al.ll201ll ). For the present system this results 



^dyn ~ 1.9 X 10" 



(6) 



which is similar to the LITE. However, as we now point 
out, a more detailed analysis shows that the ETV curve 
should be LITE-dominated. Although the harmonics of the 
fundamental frequency could arise from the eccentricity of 
one of the orbits, th ere is strong eviden ce from the radial 
velocity solution of jPerekas et al.|[201ll ) that both orbits 
are circular, which is further supported by the locations and 
shapes of the secondary minima with respect to the primary 
minima in both the close and wide orbits (see next Section). 

Accepting that both orbits are nearly (or exactly) circu- 
lar, the LITE contribution is restricted to the fundamental 
term, and there is no dynamical addition to this term. In 
this situation, the only dynamical terms that can give non- 
vanishing contributions are as follows: 



ETVdy-n 



3 mA Pi 
Stt rriAB P2 



|sin^ im sin 2 (112 — Un 



-\-- cot il sin ini [sin Uml cos 2{U2 — 1tm2) 
+ cos Zm cos Uml siu 2(u2 — Mm2)] } 



(7) 



(see Eq. (46^1 ISorkovits et all [2OO3I '). As before, indices 1 
and 2 refer to the elements of the close and wide relative 
orbits, respectively. Furthermore, im denotes the mutual in- 
clination of the two orbital planes, while Mmi and Uni2 stand 
for the angular distances of the intersection of the two or- 
bits from the plane of the sky, measured on the respective 
planes (see Fig. IA2l in Appendix We see that in the case 
of coplanarity, all these terms vanish due to sinim = 0. For 
the present situation, the second and third terms, arising 
from nodal regression (the precession of the orbital plane 
of the close pair) can also be simply omitted independently 
from the mutual inclination, due to the almost edge-on view 
of the orbital plane, as | cot ii = i cot87?6 ~ 0.02. 

As a consequence, we are in a very fortunate situation. 
Provided we accept that the 45.5-day-period sinusoidal ETV 
is caused by the above described geometrical and dynamical 
effects, the signals of the two phenomena could very easily be 
disentangled. Firstly, the amplitude of the P2-period compo- 
nent gives information about the physical dimensions of the 
close binary's orbit around the centre of mass of the triple 



^ We corrected here the erroneous negative sign in the nodal term 
(i.e. in front of i cot ii). 



system. Combining this result with radial velocity measure- 
ments of the giant companion makes it possible to determine 
the masses mA and mAB (as a function of the photometri- 
cally known sin 12), in a similar manner to a double- lined 
spectroscopic binary (SB2). Secondly, the |P2-period term 
makes it possible to determine the relative (or mutual) in- 
clination of the two orbits, i.e. the spatial configuration of 
the triplet. 

Taking into account the above considerations, the ETV 
analysis was carried out as follows. First, a general linear 
least-squares method was applied to search for the best fit 
in the following form: 

2 

f{E) ^ C0 + C1E + C2-E^ + ^ {aj sin juE + hj cosjujE) , (8) 

where the frequency was taken from the DFT analysis, and 
was held fixed. Note that its physical meaning is oj = 
where Pd and P02 stand for the eclipsing periods of the close 
and wide binaries. These quantities, strictly speaking, are 
neither equal to the anomalistic periods Pi and P2 (which 
appear in the amplitudes of the dynamical terms) [e. g. for 
7 systematic velocity Pd = Pi (l -f ^)], nor necessarily con- 
stant, especially when C2 7^ 0. Nevertheless, for our purposes, 
these differences are not significant. 

We carried out two fitting procedures: one for the com- 
plete data series, and another only for short-cadence Q7—Q9 
data. Instead of estimating and using individual measure- 
ment errors for each data points, we applied a simple weight- 
ing scheme. Namely, weights ai — OfOOOS and ai = 0"?001, 
estimated from the eclipse time determination procedure, 
were chosen for short-cadence and long-cadence minima, re- 
spectively. After a preliminary fit, points above the 3cr limit 
were removed, and the procedure was reiterated. We list our 
results from the two data sets in Table |3] while the corre- 
sponding fitted curves are shown in Fig. (2] We also show the 
phased graph in Fig. 3] The polynomial terms (i.e., ^ CiE^) 
were subtracted from this latter curve. In Table|3l along with 
the direct output of the least-squares fits, the derived physi- 
cal and geometrical quantities, and their standard errors are 
also tabulated. 

Before analysing the individual Fourier-contributions, 
we should stress, however, that there is a discrepancy of 
about 0.05 days between the wide-orbit's period obtained 
here from the LITE solution and the one determined from 
the deep eclipses directly (see later in Sect. l3.2]l . This is quite 
significant, as during the measured 17 cycle-long interval it 
would result in a shift of about 0.85 days in the occurrence of 
the eclipse events. Our light curve solution (Sect.|3| clearly 
shows that the correct period is the one obtained from the 
deep minima times in Sect. 13.21 and not the present one. 
The origin of this discrepancy is unclear. It might be caused 
by the observations of shallow minima being absent around 
the extrema of the LITE-orbit. A firm resolution will require 
further investigations on a longer time interval. Fortunately, 
this period difference is too small to infiuence the analysis 
of the Fourier terms described below. 



3.1.2 Short-period variations: light-time effect 

Considering the light-time contribution first, its most im- 
portant output is the physical size of the light-time orbit 
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we can rule out the retrograde orientation photometrically. 
Therefore, the corresponding angles are calculated only for 
prograde relative orbits. By combining the mutual inclina- 
tion, the phase term (Um2) and the visible inclination {12) 
- the latter being known from the light curve solution - we 
can calculate the complete 3D orbit of the triple system. In 
Table [3] we also give the difference of the longitudes of the 
nodes (Af2) on the sky, as well as the visible inclination i\ of 
the close system. Since i\ is also known from the light curve 
solution, this result might help to resolve the ambiguity, 
and also serves as an accuracy check for our solution. 
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Figure 4. The phased ETV curve together with the best linear 
Isq-fit for the Ql — Q9 data. Note that the quadratic term has 
been subtracted. 



of component B (at least as a function of inclination 22). 
Together with the semi-major axis of component A's orbit 
(obtained from radial velocity measurements) , this yields the 
physical masses of the wide binary (i.e., the mass of the gi- 
ant component and the total mass of the close binary). Note 
that, as one can see in Table |3] the ratio mA/rriAB has a sig- 
nificantly lower standard error than the masses individually 
and, furthermore, it does not depend on the inclination 22. 
Nevertheless, there is clearly a significant discrepancy be- 
tween the mass ratios and the masses derived from the two 
solutions. The mass ratio depends strongly on the amplitude 
of the LITE term. However, the mass of the giant component 
resulting from the pure, better-quality short-cadence data 
accurately confirms the value d erived from previous results 
and astrophysical estimations of lDerekas et al.l (|201ll ) . Con- 
sequently, in the foUowings we adopt this second (Q7 — Q9 
SC-data only) solution. 

The second parameter coming from the LITE term, the 
phase information, is less useful, but we may use it for an 
indirect checking of the accuracy of our solution. This value 
estimates a primary eclipsing mid-minimum (i.e. mid-transit 
of the small binary in front of the giant component) at BJD 
55045.21 ± 0'?16. By the use of the direct ETV-determined 
ephemeris of the wide binary (see Sect. I3.2|l we measure 
phase (f) — 0?998 for this event, i. e., the = phase oc- 
curred at BJD 55045.28762, which is clearly within the for- 
mal error. 



3.1.3 Short-period variations: dynamical effects 

Now we turn to the dynamical term. The corresponding 
Fourier coefficients (a2, 62) are almost two orders of magni- 
tude smaller than those of the LITE terms, and they are 
close to the standard errors. Consequently, the following 
results should be considered with great caution. From the 
amplitude we get sin^ im ~ 0.05, which is large enough to 
marginally verify the omission of the nodal contribution, 
but not large enough to give a numerically trustable output. 
From this result we obtain two different values for the rela- 
tive inclination. However, as will be shown in the Discussion, 



Both solutions seem to indicate a significant (13° — 15°) 
misalingnment between the two orbital planes. If this fact 
were real, a precession of the two orbital planes would oc- 
cur around the inv ariable plane of the triple system. It ca n 
be shown (see e. g. lSoderhieimll 19751 : feorkovits et al.ll2007l '). 
that the orbital inlination of the close binary would then 
vary cyclically with an amplitude of 28° — 30° on a time- 
scale of 13 — 14 years. Furthermore, the fact that the phase 
term Um2 is close to 90° or 270° (i. e. the observable in- 
clinations (ii and 12) have very similar numerical values) 
shows that this hypothetical effect would produce the fastest 
ii variations at the present epoch. This means that during 
the Ql — Q9 observational interval we should have observed 
more than 10° variation in the visible inclination (ii) of the 
close pair. This variation would have resulted in significant 
changes in the eclipse depths of the shallow minima. How- 
ever, according to our analysis (next Section) there is no sign 
of any eclipse-depth variations in the close system, and so 
we have to exclude this possibility. Consequently, the pres- 
ence of the first harmonic in the DFT-spectrum cannot be 
explained by the non-coplanarity of the orbits. 



Having ruled out both the eccentricity of the orbit (s) 
and the noncoplanarity of the orbital planes, we examined 
further possibilities by considering the effects of higher-order 
dynamical ter ms. Although al l the dynami cal te rms consid- 
ered e . g. by iBorkovits et all (|2003l . HoTJ) and lAgol et all 
(2005") disappear for coplanar and circular orbits, this hap- 
pens only within the frame of the applied approximation. 
The octuple and higher-order terms of the perturbation 
function cause non- vanishin g contributions even i n this case, 
as it was shown e. g. by ISoderhielml (|l984l '): iFord et al] 
(|2000l ). In order to check the magnitude of such forces, we 
integrated the motion numerically and calculated the sim- 
ulated times of minima. In our integration both the New- 
tonian point-mass and the non-dissipative tidal terms were 
included. The applied numeric integrator was described in 
IBorkovits et al.l (|2004l '). An analysis of the DFT spectrum of 
this higher-order, numerically-generated (and evenly sam- 
pled) ETV curve revealed the presence of the first few har- 
monics of the orbital periods at a 90% significance level. As 
the amplitudes of these peaks are lower by approximately 
two magnitudes than that of the questionable first harmonic 
in the observed curve, we can conclude that these higher- 
order effects are also insufficient to explain the structure of 
the Fourier space. Therefore, we cannot currently give any 
plausible dynamically originated explanation for the P2/2- 
period term in the ETV. 
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Table 3. Fitted and derived parameters (and their formal errors 
in the last digits) from the general linear least-squares fit to the 
ETV curve. 



Parameter 


Ql - Q9 


Q7- Q9 




0.019897(2) 


CO 


-0.000018(51) 


0.004275(1206) 


Cl 


-0.0000002(3) 


-0.0000125(40) 


C2 


0.5(4) X 10~'J 


10.8(33) X 10-8 


ai 


0.001040(28) 


0.001128(34) 


bi 


-0.001004(27) 


-0.001028(33) 


12 


-0.000001(29) 


0.000006(35) 


b2 


0.000089(26) 


0.000067(31) 


^Bab-primin [BJD] 


55051.236232(51) 


55051.240526(1206) 


Pie [day] 


0.9056768(3) 


0.9056645(40) 


APie [day/cycle] 


1.1(8) X 10"^ 


21.6(66) X 10-8 


P2e [day] 


45.518(4) 


45.517(5) 


Ob sin 12 [Rq] 


54(1) 


57(1) 


{«ab)o [°] 


-44(1) 


-42(1) 


Tab — primin [BJD] 


55045.4(1) 


55045.2(2) 


a\ sinij [Rq] 


33.43(5) 


mA/niAB 


0.617(5) 


0.629(5) 


mABsin3i2 [M0] 


4.30(15) 


4.76(20) 


niA sin^ 42 [Mq] 


2.65(10) 


3.00(13) 


niA/mAB sin^ «m 


0.042(12) 


0.031(14) 


4 [°] 


15(2) 


13(3) 


Mm2 [°] 


91(9) or 271(9) 


95(15) or 275(15) 






87.7 




88(2) or 88(2) 


87(3) or 89(3) 


An''''^ [°] 


15(2) or -15(2) 


13(3) or -13(3) 



taken from lDerekas et al.l 1I2OIII ): 

180° - Jm, 180° - ii, 180° - AO give equivalent solutions; 

fixed from the light curve solution; 

The second values arc valid for + 180°. 



3.1.4 Secular variations 

As mentioned above, the ETV curve shows weak evidence for 
continuous orbital period changes with a contant rate during 
the whole observational interval. In order to investigate this 
feature, we consider the Ql — Q9 dataset with longer time 
coverage, instead of the previously used Q7 — Q9 SC data. 
The quadratic ephemeris, calculated from this solution, for 
the shallow minima is 

MINi = 2 455 051.23623(5)+0.9056768(3)£;+0.5(4)xl0"''£;^ 

from which the rate of the constant period change is found 
to be 

^ ~ P = 2^1 ~ 0.038 s/yr. (10) 

The origin of this variation is not clear. As we men- 
tioned, any orbital precession can be ruled out due to the 
almost exact coplanarity. Due to the detached system geom- 
etry, none mass loss, mass exchange or magnetic cycles can 
be considered, as a reason. Gravitational effects induced by 
an additional, more distant and faint companion, could be 
responsible. Moreover, some interaction (e.g. tidal, magnetic 
or other) with the giant component might also be the source 
of this phenomenon. Further observations and investigations 
are needed to clarify the origin of the secular variations. 



3.2 The wide system 

For the deep minima the following linear ephemeris was 
found by a linear least-squares fit; 

MINi [BJD] = 2 455 499.9962(4) 45'?4711(2) x E. (11) 

Due to the coverage of 17 orbital cycles only, and a large 
scatter of about 0.03 days, no periodic or secular trend can 
be identified in the ETV curve. The relatively large scat- 
ter may arise from the irregular, intrinsic variations of the 
ch romospherically a ctive giant component. As it was shown 
bv lKalimeris et al.l (|2002l ). starspots can alter the measured 
mid-minimum times by ~ 0.01 days. Evidence for starspots 
(and even of eclipses of spotted regions) will be given in 
the Discussion. Therefore, we conclude that during the 2.1 
year-long observed time interval, the period of the outer or- 
bit remained constant. 



4 LIGHT CURVE ANALYSIS 
4.1 Light curve characteristics 

The fight curve of HD 181086 has at least five different com- 
ponents: 

(i)-(ii) The eclipsing features of both the close inner {Ba — Bb), 
and the wide outer (^4 — B) binary subsystems. This cate- 
gory includes not only the eclipses themselves, but also other 
effects coming from the close binarity, i. e., the ellipsoidal 
variations arising mainly from the tidally distorted shape of 
the giant component A. As we will show below, relativistic 
Doppler-beaming also produces a contribution. The reflec- 
tion effect occurs in t he close binary, bu t is negligible for 
the wide system (c. f. IZucker et al.|[2007l ). The characteris- 
tic time-scales of these variations are equal to the observed 
eclipsing periods Pi , P2 of the two subsystems. Note that the 
period ratio is almost exactly Pi ; P2 = 5 : 251, hence, in 
every fifth revolution on the wide orbit, the shallow eclipses 
occur at approximately the same orbital phases of the wide 
system. Since the shape and the duration of the deep eclipses 
are remarkably altered by the varying positions of the close 
binary members, this resonance naturally defines five differ- 
ent deep eclipse patterns (or eclipse families, which are ana- 
loguous to the Saros cycles). Furthermore, considering two 
consecutive deep primary eclipses of a given "family" (which 
occur at cycle numbers E = n and E = n-\-^, respectively), 
the intervening deep secondary eclipse of the same "family" 

(9) 

(located a.t E = n -\- 2.5) has a similar egress and ingress 
pattern, but with a 0.5 close-orbital phase shift, i. e. with 
an interchange between the shallow primary and secondary 
minima. In Fig. [5] we plotted some typical members (both 
primary and secondary) of three of the five "families" . 
(iii) The strictly periodic and regular light curve variations 
are strongly altered and distorted by irregular or semi- 
regular brightness changes with more or less similar ampli- 
tudes. This feature may come from the intrinsic variations 
of the giant primary, and suggests that this star is a chromo- 
spherically active object. Some evidence for large spots can 
be seen in the different depths and shapes of primary deep 
minima (compare Fig.[5}i and Fig. [5})): when the close binary 
transits across a darker region, the minimum is shallower. 
The irregular variation seems to be continuous, showing cer- 
tain quasi-periodicities on a 1-2 month time-scale, and could 
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Figure 5. Examples for three of the five "famihes" of the outer ecUpses. Solid curves are the raw (uncorrected) flux curves, while dotted 
ones are corrected for the intrinsic variations. Solid and dashed vertical lines denote the small primary and secondary mid-minima, 
respectively. Note the flatness of the bottom of the primary minimum-curves in panel h (especially with respect to its counterpart in 
panel a), which might be the consequence of a transit in front of a spotted region. 



have some connection with the orbital and/or rotational pe- 
riods of the giant component. 

(iv) There are further, small amplitude oscillations in the 
Ught curve with the half of the sinodic period of the close 
system with respect to the giant, which strongly indicates a 
tidal origin. 

(v) Finally, flare events were also observed during some of 
the observational runs. If these transients have their origins 
in HD 181086 then, at least in one case, we can be sure that 
it comes from the giant component, since the flare event at 
BJD 2 455 659 (in Q9) occurred during the secondary mini- 
mum of the wide system, i. e., when the close pair was totally 
occulted (Fig. (6)). 

In the present analyis, we mainly focus on the eclipsing 
features — of the light curve. As mentioned above, 
the presence of mutual eclipses in both subsystems makes 
it possible (at least theoretically) to infer some additional, 
otherwise unobtainable, physical and geometrical parame- 
ters from the hght curve solution. For example, both the fine 
structure and the variable length of the ingress and egress 
phases of the deep minima reveal information on the mutual 
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Figure 6. A possible flare event at BJD 55659 (in Q9) within 
a secondary deep minimum. The location and the amplitude of 
the eruption demonstrate clearly, that if it is a real flare event, it 
must have occurred on the giant component. 



inclination of the two subsystems in such a way that even 
the usual i, 180° — i ambiguity can be resolved, i. e. we can 
decide whether the revolutions of the two subsystems are 
prograde or retrograde relative to each other. Furthermore, 
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the combination of the shallow and deep echpses gives an in- 
dependent solution for the photometric mass-ratio in both 
the close and in the wide systems. (In Appendix [X] some 
examples are given for mining the extra information coded 
into the mutual eclipse geometry.) 

4.2 Method of the analysis 

In order to carry out this analysis, as a first step we had to 
separate the different kinds of variations in the light curve. 
While the removal of the transients (or flares) was straight- 
forward, and the small-amplitude tidally generated oscilla- 
tions do not modify significantly the eclipsing structure, the 
subtraction of the long-term intrinsic variations was a diffi- 
cult problem. We resorted to a step-by-step iterative process, 
in some steps very similar to a filtering in Fourier space. 

First, we obtained the averaged light curve of the close, 
Ba — Bb binary. Since one Kepler quarter covers ^ 100 cy- 
cles, we expect that those brightness variations which are 
independent of the close binary's orbital revolution would 
average out. We therefore binned and averaged the out- 
of-deep-eclipses parts of our light curves according to the 
eclipsing phase of the close binary. We applied this process 
for six different datasets: the three short-cadence data-series 
(Q7, Q8, Q9) were taken individually, and also together, the 
long cadence Ql — Q6 data together, and, finally, we con- 
verted the short cadence data into long cadence ones, and 
averaged the whole Ql — Q9-long LC dataset into an addi- 
tional light curve. We tried different binning numbers, and 
found 300 as an optimal solution, providing sufficient time- 
resolution and still containing enough data points in each 
cell for an effective averaging. (We have also corrected the 
phase values for LITE, although, since the cell size was ap- 
proximately equal to the full amplitude of the ETV [see the 
previous section], it had only a minor effect on the accuracy.) 
Then we obtained a light curve solution with the PHOEBE 

I I r* ( 

code (jPr sa & Zwittcr 2005) . Most of the initial parameters 
were adopted from Derckas et al.l (|201ll '). The effect of the 
giant component at this stage was considered simply (and 
crudely) as a constant third light. The initial values of this 
latter quantity were taken from the depth of the deep sec- 
ondary eclipses (where only the giant component is visible). 
In the left panel of Fig. [7] we plot the Q7 — Q9 short-cadence 
average, together with its PHOEBE solution curve. 

We also averaged the wide binary's light curve in a sim- 
ilar manner. In this case we divided one orbital revolution 
into 1000 bins (see Fig. [8]). Note that the whole Q1-Q9 time 
interval spans only ~ 17 orbital cycles, and there are also 
some gaps in the data. Therefore, we cannot expect a well- 
averaged light curve even for the full dataset. Furthermore, 
such an averaging smooths out the shoulders in the ingress 
and egress phases of the outer minima, which contain the 
most important geometric information. 

In order to recover this information, we calculated a 
preliminary net eclipsing and elliptical light curve for the 
whole triple system. For this we developed a new light curve 
synthesis code, which calculates the motions, gravitational 
interactions and mutual eclipses of the three stars simulta- 
neously. The main characteristics of our code are described 
in Appendix [B] 

For the computation of the sy nthetic curve, most o f the 
input parameters were taken from lDerekas etahl l|201ll ). re- 



fining their values with our results from the ETV analysis 
and the close binary's PHOEBE light curve solution. Af- 
ter some very minor trial-and-error fine tunings we found a 
seemingly satisfactory fit. In Fig. [8] we show two versions of 
this synthetic curve (subjected to the same averaging pro- 
cess), one including the beaming effect, and the other with- 
out. We see that the curve which includes Doppler-beaming 
(in the order of 1 ppt) gives a better fit. Despite its prelimi- 
nary stage, the fit is quite satisfactory from the first contact 
of the deep primary minimum to the next quadrature. The 
discrepancy in the other portions is probably due to the inef- 
ficiency of the averaging. An averaged residual curve is also 
shown in Fig. [H] 

As a next step, we subtracted this synthetic light curve 
solution from the raw data. This process was carried out 
individually for each quarterly dataset. The raw Q7 SC-data, 
the synthetic light curve, and the residual are plotted in the 
left panel of Fig. [1 

A discrete Fourier analysis was carried out for the resid- 
ual curves. This was applied for different datasets. First, in 
order to get the longest possible homogenous dataset, we 
made the DFT of the full Ql — Q9 LC dataset. We also made 
DFTs separately for Q1-Q6 LC data, and Q7~Q9 SC data. 
We found that the different datasets produced very simi- 
lar spectra, and consequently, similar significant frequencies. 
Using the most prominent 10-15 frequencies, we fitted sinu- 
soidal curves to the residual light curves. We found the best 
solutions, when we fitted two consecutive quarter-data to- 
gether. Finally, these Fourier polynomials were subtracted 
from the original observational data. As a final result, we ob- 
tained such a detrended 'observational' dataset, which was 
dominated by the eclipsing nature of the triple system. This 
set was used for further analysis. The step-by-step process 
for the Q7 SC-data is shown in the panels of Fig. (9] while 
three segments of the Q7 — Q9-SC DFT spectrum are plot- 
ted in Fig. 1101 The right panel of Fig. [3 showing the close 
binary's averaged light curve for the detrended Q?* — Q9* 
data, illustrates the effectiveness of this procedure. (The bot- 
tom right panel of the Figure also contains an indirect evi- 
dence for the lack of short-term variations in the inclination 
ii : a change in the eclipse depth would imply an increase of 
the point-to-point scatter during the eclipses, which is not 
seen to occur.) 

In the next stage we made a grid-search analysis with 
our code on the detrended Q7* LC-dataset. We chose this 
quarter because of its relatively regular, less-distorted shape. 
The fitted parameters were as follows: the two mass-ratios 
qi,2, the (fractional) stellar radii i?A,Ba,Bb, temperatures of 
the close binary members TBa,Bb, one of the three stellar lu- 
minosities in Kepler-band (the other two were calculated), 
the two orbital periods Pi, 2, two epochs To-1,2, two observ- 
able inclinations ii,2, and the relative longitude of the node 
of the two orbits on the sky AQ, while other parameters 
were kept as fix ones. Logarithmic limb-darkening formulae 
were applied (equivalent with Id — 2 constraint of the WD 
and PHOEBE code), with coefficients taken directly from 
PHOEBE code. The fc, internal structure constants were 
taken from the tables of lClaret fc Gimenej (|l992l ). 

In order to estimate the accuracy and reliability of the 
obtained parameters, we repeated our procedure for the 
other quarters. This enabled us to estimate the influence of 
the residual distorted, spotted features of the pre-processed 
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Figure 7. Left panel: The binned, averaged liglit curve of the Ba — Bb close binary for the Q7 — Q9 SC data (upper blue circles) with 
a typical fit yielded by PHOEBE (red line), and the standard deviations of the binned data with respect the average value of each 
individual cells (down). Right panel: The binned, averaged light curve of the Ba — Bb close binary for the detrended Q7* — Q9* SC 
data (upper blue circles) with a similarly processed typical solution curve yielded by our new synthetic code (red line), and the standard 
deviations of the binned data with respect the average value of each individual cells (down) for both the detrended observed data (blue), 
and the solution one (red). Note that the bottom curves do not represent the residuals of the upper solution curves. 
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Figure 8. The binned, averaged light curve of the AB outer 
binary for the total Ql — Q9 long cadence dataset (blue circles, 
and the synthetic eclipsing light curve averaged on the same way 
with and without Doppler-beaming (red and green, respectively). 



light curves on the solutions. All the fixed and fitted param- 
eters, as well as their estimated errors, and some derived 
quantities are listed in Table [l] 

Our final solution for Q7 data is plotted in the panels 
of Fig. [11] for some characteristic parts of the curve. 



5 DISCUSSION AND CONCLUSIONS 

We have determined a new set of physical parameters for all 
three components in the system. Our results have roughly an 



order of magnitude low er random erro r s than was achievable 
after the discovery by iDerekas et al] (|201lh . Furthermore, 
we were able to exploit the unique geometry to infer new 
parameters that were previously beyond reach. 

For the previously determined parameters, we find 
excellent agreement with the new values. For example, 
the primary's radius, combining the Hipparcos parallax 
with CHARA / PAVO onterferommetry, was measured by 
iDerekas et aD (|201lh to be Ra = 12.4 ± l.SR©. Now we 
have determined 71a = 12.46 ± O.ISR© by combining the 
stellar masses from the ETV study with the simultaneous 
light curve analysis. Similarly, the ETV analysis plus the 
SBl radial velocity measurements yielded a primary mass 
of niA ~ 3.O±O.1M0, whic h agrees with the estim ated mass 
from evolutionary tracks in IDerekas et al.l (|201lh . All in all, 
the derived physical parameters draw a consistent picture of 
the system, proving that despite the difficulties in the light 
curve modelling, our method yields robust results. 

A preliminary compa rison with models fro m the BASTI 
JPietrinferni et al.] |2004| ) and Dartmouth (jDotter et al.l 
l2008h databases shows that the fundamental properties for 
all three components are consistent with solar-metallicity 
isochrones with ages ~ 300 — 500 Myr, although the dwarf 
radii appear to be significantly larger than expected. More 
detailed comparison using the near model- independent prop- 
erties presented here will allow powerful tests of stellar evo- 
lutionary theory, such as tidal effects on the mass-radius 
relati on for low-mass s tars in close-in binary systems (see, 
e. g.. lKraus etHlbOllI ). 

One important question in relation to the giant pri- 
mary is its evolutionary stage, being located in a part of 
the H-R diagram where JZ-shell-burning stars ascending the 
first red giant branch overlap closely with _ffe-core-burning 
giants (in other words, there is an age uncertainty that can- 
not be resolved from the evolutionary tracks alone). Dy- 



© 2012 RAS, MNRAS 000.[llf23l 



16 



T. Borkovits et al. 



1.010 
1.005 

; 1.000 

I 0.995 

1 0.990 

: 0.985 

0.980 

0.975 
; 0.010 
: 0.005 
1 0.000 
i -0.005 
-0.010 




H0 '*«.» 



1.010 

1.005 

J 1.000 

"g 0.995 

1 0.990 

^ 0.985 

0.980 

0.975 
i 1.010 
"S 1.005 
% 1.000 
S 0.995 
^ 0.990 





1 Ak 


<j ' I ■ 



0.990 

; 0.010 - 
: 0.005 

I 0.000 
i -0.005 
■ -0.010 



55460 55480 55500 55520 55540 55560 
BJD 



55460 55480 55500 55520 55540 
BJD 



55460 55480 55500 55520 55540 55560 
BJD 



Figure 9. The process of the removal of the intrinsic Hght curve variations from the raw data for the Q7 SC observations. Left panel: 
The subtraction of a preliminary synthesized eclipsing light curve (green) from the original Q7 data (upper red curve) results a residual 
curve of the irregular variations (lower red curve). Middle panel: After a DFT-search of the significant frequencies in the residual curves, 
the intrinsic variations are represented by the corresponding Fourier polynomial (green), and this latter curve was subtracted from the 
original data (upper red). The detrendcd Q7* data are plotted in the middle lower panel with red color. Right panel: The final light 
curve solution (green) was fitted to this Q7* dataset (upper red). The residual curve can be seen in the bottom panel. 
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Figure 10. Three different frequency-regions of the DFT spectra obtained for three different datasets of Q7 — Q9 SC data. Red curves 
denote the spectrum of the raw data, green ones represent that of the synthetic light curve solution, while blue spectra refer to the 
residual light curve. Dashed lines in the left panel show the frequencies belonging to the eclipsing period, and its first harmonic. Note 
the different scales along the Y-axes. See text for details. 
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Figure 11. Three zoom-ins into the Q7 and Q8 light curves with solutions, and residuals. Upper panels: Points show the original Q7 and 
Q8 data, while dashed lines represent the sum of the solution for the detrended Q7* and Q8* data, and the Fourier-modelled intrinsic 
variations. Bottom panels: The residual curves. The left and middle panels show small parts of the Q7 data just after a large primary 
minimum, and in the following quadrature phase, respectively, while the right panel represents the surroundings of a large secondary 
minimum from the more distorted quarter Q8. 
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Table 4. Stellar and orbital parameters derived from the com- 
bined ETV and synthetic light curve analysis. (The numbers in 
parantheses are the estimated errors in the last digits.) 

orbital parameters 

subsystem 





Ba-Bb 


A-B 


P[d] 


0.9056768(2) 


45.4711(2) 


Tmini [BJD] 


2455051.23623(5) 


2455499.9962(4) 


a [Kq] 


4.777(39) 


90.31(72) 


e 


0.0 


0.0 








i [deg] 


86.7(14) 


87.5(2) 


An [deg] 


0.0(5) 




«m [deg] 


0.8(14) 




1 


0.95(3) 


0.595(5) 


Lsec/ L'YOT 


0.3468 


0.0078 


stellar parameters 




Ba Bb 


A 


fitted and/or derived parameters 


relative quantities 


^pole 


0.1798 0.1664 


0.1376 


^side 


0.1808 0.1672 


0.1379 


^point 


0.1826 0.1687 


0.1382 


^back 


0.1822 0.1684 


0.1382 


absolute quantities 


m [Mq] 


0.915(34) 0.870(43) 


3.0(1) 


R [Rq] 


0.865(10) 0.800(20) 


12.46(15) 


Teff [K] 


5100(100) 4675(100) 


5100(100) 


Lhoi [Lq] 


0.447(37) 0.270(27) 


92.812(7615) 


logg [dex] 


4.53 4.58 


2.73 


fixed quantities 


k2 


0.020 0.020 


0.033 


/3 


0.32 0.32 


0.32 


A 


0.5 0.5 


0.5 


Xhol 


0.71476 0.71476 


0.71159 


2/bol 


0.13026 0.13026 


0.12561 


XK 


0.70835 0.70835 


0.70074 


VK 


0.16354 0.16354 


0.16609 



namical considerations can help here, too, via comparing 
the orbital configurations with theoreti cal tidal circulariza- 
tion ti me-scales. According to Eq. (7) in lVerbunt fc PhinnevI 
l|l995h . which was based on the works of lZahnI l|l977l . ll989l ). 
a binary with the same parameters as HD 181068 A and 
B (=Ba-|-Bb) is expected to be circularized under a period 
limit of Pcirc ~ 15 days for //-shell burning primary. With 
the observed P2 ^ 45 days and the perfectly circular orbit, 
theory implies indirectly that the primary must be older, so 
that in the //e-core burning phase. The question, however, 
is more complicated because of the binary nature of the sec- 
ondary. This causes additional complications by the tidal 
oscillations that are expected to affect the convective enve- 
lope of the primary. It is not known if the tidal damping is 
effective enough to shorten significantly the circularization 
time. 

Considering the other orbital parameters, our solution 



for the orbital inclination of the close binary (ii = 86?7±1?4) 
has a relatively large uncertainty. This is not surprising be- 
cause, being a partially eclipsing pair, the observable incli- 
nation is very sensitive to any additional third hght, i. e. in 
the present case, for the hght of the giant primary, and par- 
ticularly, for the continuous variation of this extra amount 
of light. The inclination of the wide system was found to 
be (22 — 87°. 5 ± 0?2). From these two values alone, in the 
absence of any other information, we would be able to say 
nothing about the mutual inclination of the system. How- 
ever, the simultaneous hght curve fit of this triple eclipsing 
system provides a direct and powerful method for determin- 
ing this quantity. As mentioned above, this comes from both 
the fine structure and the timings of the ingress and egress 
phases of the deep eclipses. This is illustrated in Fig. 1121 
where we plotted the first two deep eclipses of the Q7 SC 
light curve. The only difference between the different col- 
ored curves is the AQ parameter, and consequently, the 
mutual inclination. While the out-of-deep-eclipse parts, and 
the totality-of-eclipse periods of the fight curves are iden- 
tical, the ingress/egress fine structures, and the moments 
of the contacts differ significantly, and this makes the Afl 
adjustable parameter (and so the mutual inclination) to a 
well-determined quantity. Furthermore, even the im = 1° 
curve is definitely separable from its retrograde counterpart 
im = 179° . (Note, however, that this separation is only pos- 
sible when the masses or the radii are different in the close 
binary.) A combination of the obtained Afl parameter with 
the two observable inclinations results in a mutual inclina- 
tion of im = 0?8 ± 1?4, which suggests an exact coplanarity. 
This is in accordance with the lack of the eclipse-depth vari- 
ation of the shallow eclipses. 

Finally, we briefly comment on the other features of the 
light curves. First we consider the irregular, or semiregu- 
lar brightness changes, which likely originate from chromo- 
spheric activity. Evidence of the presence of spotted regions 
on the giant's surface was shown in the previous section 
(see e. g. Fig. Further characteristics can be deduced 
from the comparative investigation of the DFT spectra of 
the raw observed hght curve, the synthetic and the residual 
ones (Fig. [To}. What can be seen well even at the first glance 
is that in the low frequency domain (left panel), the spec- 
trum of the observed data remarkably departs from that 
of the synthetic data. While in the synthetic eclipsing, el- 
lipsoidal data the dominant frequency corresponds to the 
half eclipsing period of the wide system, the highest peak 
of the original data is located about the eclipsing period it- 
self. Furthermore, this latter peak is clearly a double one, 
whose two peaks are already well separated in the spec- 
trum of the residual light curve (i. e. after the removal of 
the eclipsing and ellipsoidal features). In our interpretation 
these two peaks might have a rotational origin. The good 
correspondance of this pair of peaks with the orbital pe- 
riod proves the synchronised rotation of the primary, while 
its splittin g might giv e an evidence of differential rotation 
(see e. g. [oiah et a l.i |200d). Note, that the spectrosc opi- 
cally obtained Vrot sini — 14 kms~^ (jPerekas et al.|[201ll ) for 
Ra = 12.5 Rq, and sini2 = 87°. 4 result in Prot = 45'? 474, 
which is also in very strong correspondance with this result. 

Considering the high-frequency end of the DFT spec- 
tra (right panel of Fig. I10|l . three distinct peaks can be 
identified at /i = 2.20829 d'^, /2 = 2.16431 d"^ and 
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Figure 12. Synthetic light curve for an outer secondary (left), and a primary (right) minimum, calculated with different mutual 
inclinations. See text for details. 



/s = 2.12032 d~^, from which three, the first is exactly 
the half of the eclipsing period of the close binary, while 
the other two are /2 = /i — /o, and /s = /i — 2/o, where 
/o = 0.04398 d~^ corresponds to the half of the eclipsing 
period of the wide system. Such a way, the tidal origin of 
this small amplitude oscillation on the surface of the giant 
primary is out of question. These oscillatory features will be 
investigated in details in a forthcoming paper. 
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APPENDIX A: DETERMINATION OF SYSTEM 
PARAMETERS FROM THE MUTUAL 
ECLIPSES 

In this Appendix we show examples of how several system 
parameters can be determined from the geometry of the 
large, mutual eclipses. Strictly speaking, the most impor- 
tant condition for the validity of the following calculations 
is not the mutuality of the shallow and the deep eclipses 
themselves, but rather the fact that due to the hierarchical 
configuration of the triple system, the deep eclipses contain 
some mixtures of individual eclipses of the two members of 
the close pair with respect to the more distant giant com- 
ponent, which produce small changes in the deep eclipse 
configurations from eclipse to eclipse, and even between the 
ingress and egress phases of the same event. Our algorithm 
is a natural extension of the well-known methodology of de- 
termination of the relative radii of the stars (with respect to 
their separation, and as a function of their orbital inclina- 
tion) purely from the eclipse geometry, commonly used from 
the very beginning of eclipsing binary studies. 

The usual method in binaries with spherical compo- 
nents is well known: the sky-oriented distance of the stellar 
disc centers is R\ + R2 at the first and last contacts (i. e., 
at the start of ingress and at the egress phases), and if the 
eclipses are total (either transit or occultation), the same 
distance is R\ — R2 at the second and third contacts (i. e., at 
the end of the ingress and at the start of the egress phases). 
Then expressing the projected distances with the orbital el- 
ements and time, and measuring the eclipse durations (both 
the one from the first to the last contacts, and the totality 
length from the second to the third contacts), the individual 
fractional radii of the stars can be determined. 

For our triple star configuration, the egress and ingress 
phases of the deep eclipses show a complex pattern. The 
two dwarf members of the close binary may enter in front 
of or behind the giant's disk individually, or even simul- 
taneously (see Fig. lAlfl . Additionally, during an entry the 
stars' velocities, directions and distances (both physical and 
projected) relative to the giant component change continu- 
ously, producing variable length and shape in the egress and 
ingress patterns. Anyhow, no matter how complex an egress 



or ingress pattern is in itself, every eclipse event contains one 
and only one first, second, third, and fourth contacts. And 
furthermore, assuming that a given contact is not strongly 
altered by a just ongoing shallow eclipse event, we can sim- 
ply and unambiguously decide which member of the close 
binary takes part in the given contact event. For example, 
in case of prograde revolution, the very first contact of a pri- 
mary transit is produced by the eclipser of the last shallow 
eclipse event, i. e., if the last event was a small secondary 
minimum, then the very first contact of the large primary 
transit is produced by the primary of the close pair. 

Let us consider the projected distances at the disk cen- 
tres in the moments of the contacts. In the present situation 
the projected distance between the eclipser and the eclipsed 
stars no longer will be the projected radius vector of a Ke- 
plerian relative orbit, but comes from the superposition of 
two Keplerian orbits: the absolute orbit of the close binary 
members around their center of mass (CM), and the rela- 
tive orbit of this CM around the giant component. The most 
convenient and practical description of the present scenario 
uses Jacobian vectors. The first Jacobian vector (pi) is di- 
rected from rriBa to meb, i- e., it is the radius vector of the 
close binary's relative orbit, while the second one (p2) orig- 
inates from the CM of the close pair, and ends in ttia, i. e. 
it is the radius- vector in the wide pair (see Fig. IA2|I . With 
these notations, the position vectors connecting the three 
stars mutually are 

dBaBb = Pi, (Al) 
dsaA = P2 + -- ^Pl, (A2) 

l + qi 

ffebA = p2 - — Pi, (A3) 
1 + <?i 

where, as before, qi denotes the mass ratio of the close pair. 
In the usual astrometric frame of reference the right-handed 
X and y coordinate axes lie in the plane of the sky, while the 
z axis points outwards from the observer. In the astrometric 
convention, x points to the celestial north pole. In case of 
photometry, however, both the eclipsing light curve, and the 
radial velocity are invariant with respect to any rotation in 
the plane of the sky, and so, in the absence of any additional 
information on the spatial orientation of the intersection of 
the orbital plane and the sky (i. e., fi), we are free to use 
any orientation for the x axis. In the context of modelling 
eclipsing binaries, the coordinate equations take their sim- 
plest form if one of the axes in the plane of the sky coincides 
with the nodal line. In this case, the other axis gives the 
direction of the projected orbital angular momentum vec- 
tor. In the present case, however, we cannot use this latter 
formal simplicity, because of the differing orbital planes of 
the close and the wide orbits. 

It is well known from the textbooks of celestial mechan- 
ics and/or astrometry that in such a frame of reference the 
cartesian coordinates of a Keplerian orbit can be written as 

X = r-[cos(« -f oj) cos f2 — sin(?; + cj) sin f2 cos j], (A4) 
y = r-[cos(« + cj) sin f2 + sin(u + oj) cos f2 cos j], (A5) 
z = r sm{v + uj) sini, (A6) 

What is important for us is the projected distances onto 
the plane of the sky, instead of the spatial ones. In vectorial 
forms e. g. 
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Figure Al. The "walzer" of the close pair in front of (left), and behind (right) the giant primary, projected on the sky. The upper panels 
are identical with Figs. [5j and a, respectively. The dashed horizontal lines denote the Rp^ ± -Rea.b distances from the center of mass of 
the giant, i. e. the outer and inner contact places. Vertical lines connect the moments of the different contacts with the corresponding 
light curve points. 
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1 + gi P2 
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(^1 — sin'^ i\ sva^ ui) 



,(A10) 



where Ui = Vi + uji gives the true longitude of the given 
object measured from the node and furthermore. 



A = cos wi cos W2 + sin mi sin W2 cos in 



(All) 



is the direction cosine between vectors pi and p2, in which 
expression wt = m — Umi denotes the true longitude mea- 
sured from the intersection of the two orbital planes, while 
itnii is a nodal longitude-like quantity, namely the angular 
distance of the intersection of the given orbital plane from 
the sky (cf. Fig. IA2|) . It can also be seen in this figure that 
the three inclinations form angles of that spherical triangle, 
the sides of which are the three node-like arcs Af2, itmi and 
Uin2- Consequently, the two observable inclinations (ii, 12) 



and the difference of the nodes of the close and wide orbits 
(AJl — ^l2 — i^i) unambiguously determine the remaining 
quantities (i. e., im and Um-s) with the copious identifica- 
tions of the spherical triangles, from which some of the most 
useful ones in the present context are as follows: 

cos in-i ~ cos ii cos 12 -f- sin ii sin 12 cos AQ, (A12) 

sin im cos itni2 ~ — cos ii sini2 -I- sin ii cos 42 cos An(A13) 

sininisinum2 ~ siniisinAf2, (A14) 

cosii — cosi2 cosim — sini2 sinini cositm2,(A15) 

cos Uml ~ cos Um2 COS AQ + sin Um2 sin Af2 cos 12 . 

(A16) 

At this point we note that in case of coplanarity 
Eqs. (|A9IA10|) become more simple, not only due to sin ii = 
sin 12 , but also because then the direction cosine A is simply 



A = cos(u2 ~ Ul). 



(A17) 



In what follows, we assume that both the inner and the 
outer orbits are circular, as it happens to be in HD 181068. 
In this case pi ^2 = a 1,2, which means a substantial simplifi- 
cation in our treatment. At this point we are in the position 
to give the functional dependencies of the stellar sizes from 
the orbital elements. Namely, for any contact of the deep 
eclipses 

Ra ± RBa,Bb = a2fBa,Bb{ii,i2,AQ,qi,ai/a2;ui,U2), (A18) 

where the plus and minus signs hold for outer and inner con- 
tacts, respectively. Moreover, for the partial shallow eclipses 
for the outer contacts we can also write 

Rbu + RBb ^ aifBab{ii;ui). (A19) 

In these equations the independent (time-like) variables are 
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away from the observer 




n,-n, 



Figure A2. An illustration of the angles and other quantities 
used in the text. 



hidden in the Ui-s, while the other parameters are constants. 
The Ui longitudes are very close ly related to the eclipsing 
phase s. It is well known (see e. g. iGimenez fc Garcia-Pelavol 
Il983l ) that for circular orbits in the moment of a mid- 
minimum u = 90° or It = 270°. Considering the shallow 
minima, in the present situation, as pi points toward the 
secondary of the close pair, i. e., ui refers to the relative or- 
bit of the secondary, or _B6-component, therefore ui — 90° 
for the secondary minimum, and ui = 270° for the primary 
one. So the connection of ui with the eclipsing phase of the 
close pair is very simple, ui = 360° x (pi —90°, or converting 
it to time directly, ui{t) — j^{t — To) — ^ , where To is a 
mid-primary minimum time. 

The calculation of it2 requires a bit of extra care. First, 
we have to keep in mind that p2 is oriented from component 
B (i. e. the centre of mass of the close pair) towards com- 
ponent A, and so in the present formulae U2 refers to the 
relative orbit of the main, giant component A around the 
smaller and fainter component B, i. e., the secondary of the 
wide system. Consequently, in this case U2 — 90° formally 
means not the secondary but the primary deep minima. Fur- 
thermore, the determination of a reference mid-minimum 
time is not so simple. Because the positions of the close bi- 
nary members on their orbits are different at the beginning 
and at the end of a deep eclipse, the mid-minimum does 
not occur exactly at half-time between the first and the last 
contacts (similarly to the eccentric case). Nevertheless, by 
the use of an averaged fight curve (like the one in Fig. [S| , or 
of a radial velocity curve, or from the average of several ap- 
proximately determined mid-minima times we can obtain a 
satisfactory reference mid-minimum moment, and then U2{t) 
can be calculated in the same way as ui. 

In the next step we use the fundamental difference be- 
tween a traditional simple eclipsing binary and our triple 
system. In a single eclipsing binary all the eclipses are sim- 
ilar, i. e., all contacts occur at the same orbital phase, lon- 
gitude (u) and, furthermore, for circular case the eclipses 
are geometrically symmetric in time around their midpoint. 
Due to the latter, for the first and fourth (last) contacts 
sin'^ ui = sin^ uiy , and a similar relation can be written for 



the inner contacts. As a consequence, we have only one equa- 
tion for ^^'^^^ and one for El^zHs.^ and so, some extra way 
is needed to resolve the inclination dependence. In opposi- 
tion, in case of our deep eclipses, the configuration of the 
first and last contacts, and the second and third contacts 
as well, generally vary from eclipse to eclipse, and they are 
even diflerent within the same event. Consequently, we can 
get separate sets of equations ()A18|I for different ui and U2, 
from which the unknown parameters can in principle be de- 
termined by numerical methods. 

Furthermore, once ai/a2 is known, q2 can also be cal- 
culated easily by the use of Kepler III: 

ITS = 

The individual relative radii of the three stars can also be 
derived from the combination of Eqs. (|A18p . even without 
the use of Eq. (|A19|) . On the other hand, these individaul 
radii can be determined also in the case when both the deep 
and the shallow minima are partial. 

Several difficulties arise, however, during the practi- 
cal application of this method. For example, in our case of 
HD 181068, the light curve distortions of non-eclipsing origin 
cause difficulties in the accurate determination of locations 
(and so times) of the contacts. Furthermore, our stars are 
not exactly spherical, and finally, due to the 5 : 251 mean- 
motion resonance, we can get only a limited number of dif- 
ferent eclipse configurations. On the other hand, there are 
some additional results that may serve as auxiliary sources 
of information. For example, Eq. (|A19|) could be used as 
an additional equation for ii , or even ai /a2 . Radial velocity 
and/or ETV results provide further constraints or equations. 

In the foUowings we give a practical example for 
HD 181068. To do this, we use our simulated fight curve solu- 
tion. By this trick we avoid the practical problem of correct 
and accurate identification of contact times, since our pur- 
pose is simply to demonstrate the theoretical eflectiveness 
of this method. Furthermore, our finding about the orbital 
coplanarity makes our formulae as simple as possible, there- 
fore the whole calculation can be performed analytically. 

In Table lAll we list the available contact moments for 
the large secondary eclipses occuring within the Q7 — Q9 
SC data. We used only secondary occultations to ignore the 
additional uncertainty arising from the limb-darkening dur- 
ing the primary transits. According to the last column (i. e., 
U2), one can see that its value may vary by a few tenths of 
degree for both the inner and outer contacts. Although this 
may look like a small variation, note, however, that a typical 
shift of 0?5 in U2 translates to a change in the occurrence of 
the corresponding event of St ~ 1.5 hours. 

According to the above table, we have 3 different outer 
and 3 different inner contact moments for star Ba, while for 
Bb these numbers are 2 and 5, respectively. We say different, 
as events E — —0.5 and E = 4.5 belong to the same eclipse- 
family, and consequently, the corresponding moments are 
so similar to each other that we counted them only once. 
However, even in this case, we can write more equations 
than what is necessary. 

As an example, we show a concrete calculation. Using 
the first and last contacts of the E = —0.5 and the first 
contact of E — 0.5 events, substituting the corresponding u 
values into Eq. (|A9I) . and subtrating their squares from each 
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Table Al. Contact times for outer secondary eclipse events 



Table A2. Analytic results from eclipse geometry and the origi- 
nal values 



No 


COHtclCt 




MBJD 




112 


-0.5 


I 


Ba 


55476.1096 


313? 75441 


260? 89609 




II 


Bb 


55476.4245 


78? 92486 


263? 38919 




III 


Ba 


=;^477 Qfi77 


332°33560 






IV 


Ba 


'^'^A7H 4799 






0.5 


I 


Ba 


55521.5177 


3? 14288 


260?39810 




JJ 


Ba 






9fi4°4'^7A9 


1.5 


III 


Bb 


55568.9434 


134? 51263 


275?87372 


2.5 


I 


Bb 


55612.4733 


157? 33030 


260? 50577 




II 


Bb 


55612.9903 


2? 83403 


264? 59893 




III 


Bb 


55614.3571 


186? 12706 


275? 42007 


3.5 


II 


Bb 


55658.3965 


51?46726 


264? 08590 




III 


Ba 


55659.9241 


298? 67710 


276? 18012 




IV 


Bb 


55660.2422 


65? 11954 


278? 69856 


4.5 


I 


Ba 


55703.4516 


320? 54080 


260? 79316 




II 


Bb 


55703.7629 


84? 28028 


263? 25777 




III 


Ba 


55705.3125 


340? 23499 


275? 52616 




IV 


Ba 


55705.8234 


183?31398 


279? 57102 



other we get three different equations, from which we need 
only two. As such, 



sin^ i (q2 + /32a; + 723;^) + ^22; = 0, 
sin^ i (^a-j + I3:ix + jsx^^ + S-iX = 0, 



where 



X — 

and 

Q!2,3 
/92,3 
72,3 
S2,3 



ffl2 



(sin^ U2)i - (sin^ 1^2)2,3, 
2[(sinui)i(sinu2)i - (sin ui)i(sin 1*2)2,3] 
(sin^ ui)i - (sin^ 1^1)2,3, 
2[cos(jii — 112)2,3 — cos(iii — 112)1, 



(A21) 
(A22) 

(A23) 



(A24) 
(A25) 
(A26) 
(A27) 



respectively. Eliminating the (5-terms and using the fact that 
sin^ i 7^ 0, we get a second order equation in x, i. e. the 
orbital ratio. Obtaining x, the inclination can simply be cal- 
culated from any of the two equations. A similar treatment 
can be applied to the other component Bb. In this case, as 
we had only two outer contact times, we used the first two 
inner contact moments for the second equation. (Note that, 
for Bb, the signs of /? and 5 coefficients should be changed!) 
When the two ratios aBa/a2 and a-BhO.2 are known, both 
qi and ai/a^ can be immediately calculated, and then 52 
follows too. Finally, the fractional radii of the three stars 
are also easily detemined. Remaining at the present sample, 
we first determine Ra and 7?Bb from the two (one outer and 
one inner) contact equations of component Bb, and then i?Ba 
can be calculated from any of the outer contact equations for 
Ba-star, without using any inner contact moments. (On the 
other hand, we can also calculate inner contact moments for 
component Ba, of course.) There is, however, another pos- 
sibility to determine all these quantities without using inner 
contact times. This is because we can also use the outer 
contact-equation of the shallow eclipses. Thus the theoret- 
ically minimal data needed are the moments of: one outer 



parameter 


contacts used 




calculated 


original 


"Ba 
"2 


(I,IV)_0.5,l0.5 




0.023 


0.026 


i 


(I,IV)-0.5,l0.5 




83? 7 
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°Bb 


II_0.5,IIIl.5,l2. 


,5,IV3.5 


0.027 


0.029 


i 


II_0.5,IIIl.5,l2. 


,5,IV3.5 


87° 


86? 7/87? 5 


91 






0.85 


0.95 








0.049 


0.053 


92 






0.439 


0.545 


RA/a.2 


from Ba 




0.170 


0.138 


RA/a.2 


from Bb 




0.139 


0.138 


RBa/a.2 






0.009 


0.010 


RBb/a2 






0.010 


0.009 



contact for a shallow eclipse, three outer contacts for deep 
eclipses of one of the components, and two outer contacts for 
the other component. This means that we do not need any 
inner contact times (so the method works for partial eclipses, 
too), the usually better measurable outer contact times are 
sufficient. In Table \K2\ we give our results. For comparison 
we also give the corresponding parameters of the synthetic 
light curve. 

Finally, combining these results with the LITE solution, 
which returns ob sini in physical units (see Sect. 2), we can 
also calculate all the masses and stellar and orbital sizes in 
physical dimensions. So, we can conclude that in the case of 
triply eclipsing hierarchical systems (i. e., where all the three 
objects eclipse each other at least partially, but not necessar- 
ily simultaneously) a high-precision single-band photometry 
of the eclipses is at least in principle efficient for determining 
all the above presented quantities in physical units. 

(The above calculations were made for coplanar and 
circular orbits. In the non-coplanar case our equations can 
be numerically solved in a similar way. The eccentric case is 
more complicated, but the asymmetry of the eclipses both 
in length and in phase gives all the required information too, 
so the difficulty is only practical.) 



APPENDIX B: LIGHT CURVE SYNTHESIS 
CODE FOR HIERARCHICAL TRIPLE 
SYSTEMS 

The code is largely based on the well-known Wilson- 
Devinney program, whi ch is being continuously de veloped 
version 



from its first version (Wilson & Devinnev 'l97l|) up to 
now (Wilson 2008; Wilson & Van Hammc 2009). (See also 
iKaUrath fc Milond 120091 . C hapters V I and VII.). The 
PHOEBE Scientific Reference (|Prsall2006l ) was also used as a 
cook book. Some of the subroutines were borrowed directly 
from the Fortran code of the WD program (converting them 
from Fortran to C). However, a number of significant alter- 
ations were also applied. First, our code calculates the mo- 
tion and positions of all the three stars, and naturally, the 
mutual eclipses (i.e. when an eclipse event of the close bi- 
nary occurs in front of the disk of component A, or during 
the egress or ingress phase of the wide eclipse events). The 
mutual tidal interaction of the three stars is computed for 
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every moment. In order to do this in a more simple way, 
instead of the usual two-mass point Roche-model, the stel- 
lar surfaces (and the local gravities as well) were calculated 
from the first order (linear) series expansions of the potential 
of a moderately distorted spherical body. In this formalism 
the stellar radius can be written in the following form: 

= i? 1^1 + 5]/^ ' (Bl) 

where the amplitudes of the first order tidal distortions 
caused by star k on star i are 

while the amplitude of the rotational distortion of star i is 

In the equations above Ri stands for the undistorted radius, 
kj denotes the j'-th internal structure constant of star i, 
pik the distance of the two stars, LJi the rotational angular 
velocity of the star, while the direction cosines in the argu- 
ments of the given Legendre polynomials Pj are the angle 
between the radius vector of the given surface element and 
the axis of the tidal bulge (practically the radius vector con- 
nects the centre of mass of the two stars) in the tidal terms 
(A"), and the angle between the s ame surface element and 
the axis of stellar rotation iy'). See lKopall (|l97^ . Chapter II 
for details. 

Strictly speaking, for strongly distorted systems this 
formalism is less accurate than the closed form of the Roche- 
model, but in the present situation, due to the moderate 
distortion of the present stars, it is adequate. Furthermore, 
besides the obvious advantage coming from linearity, it also 
treats the stars more realistically, as it no longer attributes 
infinite central mass densities to them. 

Another improvement is the inclusion of the relativistic 
Doppler-beaming effect. The contribution to the total beam- 
ing is calculated for each surface cell of the three stars in- 
dividually, and the radial velocity contribution coming from 
the stellar rotation is also taken into account, so theoreti- 
cally the code is able to model the beaming analogous of the 
Rossiter-McLaughlin effect too. (Regardless, in the present 
system this is insignificant.) The beaming effect is well illus- 
trated in Fig. [HI 

Finally, we have also included the light-time effect. It 
is applied only to the wide subsystem, meaning that the 
positions of component A and component B were calculated 
in different moments according to their different distances 
from the observer, and furthermore, a second time-delay was 
also calculated in modelling the momentary tidal fields effect 
on each component. 

In its present version, the program has 3 x 11 star- 
specific global physical parameters (masses, radii, tidal dis- 
tortion [k2.A\ parameters, effective temperatures, chemical 
abundances, gravity darkening, two bolometric limb dark- 
ening coefficients, and bolometric albedos), 3x3-1-1 filter- 
specific parameters (luminosities, two limb darkening coef- 
ficients, and fourth light), 2x6-1-1 orbital parameters (the 
six orbital elements of the two orbits and the systemic ra- 
dial velocity), and finally, 3x6 Eulerian angles and angular 



velocities describing stellar rotations. There are also several 
flags which turn on and off various constraints between dif- 
ferent variables. Some of them are identical with those used 
in the WD program, but there are additional ones, due to 
the specific model. For example, instead of masses (which 
are undefined in the case of a simple two-body eclipsing 
light curve), one mass (usually ttia), and two mass ratios 
{q\ and 92) can also be used as input parameters. Similarly, 
instead of absolute radii, the use of fractional radii is more 
practical for light curve solution, although if the mass of the 
tertiary (rriA) and the outer mass-ratio (§2) is known and 
fixed from ETV-solution, and of course, the orbital periods 
are also fixed, the use of absolute or fractional radii is fully 
equivalent. 
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